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Abstract 

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of 
genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via 
selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 
Chinese Holstein cows with the PLINK software. Genotypes were obtained with the lllumina BovineSNPSO Bead chip and a 
total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide 
significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome 
regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of 
these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant 
associated with multiple SNPs, especially CI 8:0 (105 SNPs), CI 8 index (93 SNPs), and CI 4 index (84 SNPs); Several SNPs are 
close to or within the DGATl, SCDl and FASN genes which are well-known to affect milk composition traits of dairy cattle. 
Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising 
candidates for C10:0, C12:0, C14:0, C14:l, C14 index, C18:0, C18:ln9c, C18 index, SFA, UFA and SFA/UFA were found, which 
composed of HTRIB, CPM, PRKGl, MINPP1, LIPJ, LIPK, EHHADH, MOGATI, ECHSl, STATl, SORBSl, NFKB2, AGPAT3, CHUK, 
0SBPL8, PRLR, IGFIR, ACSL3, GHR and OXCTl. Our findings provide a groundwork for unraveling the key genes and causal 
mutations affecting milk fatty acid traits in dairy cattle. 
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introduction 

Fat is the major energy substance in milk and more than 50% 
milk total energy comes from milk fat, which accounts 3-5% of 
mUk contents. Fat nutrition value depends on fatty acids. 
Monounsaturated fatty acids (MUFA) have a favourable effect 
on human health because of its cholesterol-declining properties. 
Polyunsaturated fatty acids (PUFA) of the n-6 and n-3 series are 
essential nutrients that exert an important influence on plasma 
lipids and serve cardiac and endothelial functions for prevention 
and treatment of coronary heart diseases [1]. Conjugated linoleic 
acid (CLA) has effects on bone formation and the immune system 
as well as fatty acids and lipid metabolism and gene expression in 
numerous tissues [2]. Saturated fatty acids (SFA) lead to increase 
the concentration of low density lipoprotein (LDL) cholesterol and 
cause cardio cerebral vascular disease [3]. Therefore, changing the 
proportions of dietary fat by decreasing SFA and increasing 
MUFA and PUFA is vital to Human health. It is suggested that the 
ideal balance would seem to approximate 1:1.3:1 for SFA:MU- 
FA:PUFA [4]. 



From the genetics point of view, milk fatty acids are complex 
traits influenced by non-genetic factors, such as breed, herd, stage 
of lactation, etc [5,6] and genetic factors [7]. Bovine mUk fatty 
acids have been found to be heritable, with heritabUity estimates 
ranging from 0.22 to 0.71 [8,9]. Short and medium chain C4 to 
C16 saturated and monounsaturated fatty acids, which are 
synthesized de novo in the mammary gland, have moderate to 
high heritabUity (0.4-0.6) [8,9] . Long chain fatty acids (above C 1 6) 
are derived from circulating plasma lipids, whereas have low to 
moderate heritability (about 0.2) [8,9]. Identifying genes and loci 
responsible for the genetic variation is expected to contribute 
greatly to our understanding of milk fatty acids synthesis, and to 
develop a marker-assisted selection to improve fatty acids in dairy 
breeding program in future. In the past few years, candidate gene 
and quantitative trait locus (QTL) mapping approaches have been 
implemented to detect genes or QTLs for milk fatty acid traits. A 
few promising loci, e.g. DGATl p.Lys232Ala and SCDl p.Ala293- 
Val [10,11,12] and a large number of significant or suggestive 
genomic regions [13,14] were identified. Although the above two 
methods have got a few prominent findings, identification of causal 
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Table 2. Genome-wise and suggestive significant SNPs for short- and medium-chain saturated fatty acid traits (SCFA and MCFA). 



Trait 


Rank' 


SNP name 


Chr. 


Position(bp) 


l^cdicsL Ucilc/ 

Candidate Gene 


(bp) 


Raw P_value 


Bonferroni 


CI 2:0 


17 


BTA-1 20369-no-r5 


0 


0 


NA 


NA 


1 .47E-05 


0.599 


CI 4:0 


19 


BTB-00046603 




103185339 


SI 


Within 


1 .06E-05 


0.432 


CI 4:0 


26 


BTB-01201574 




109284636 


SCHIPl 


Within 


1 .84E-05 


0.745 


CI 0:0 


4 


ARS-BFGL-NGS-23583 




132049236 


A4GNT 


14386 


4.10E-06 


0.166 


CI 2:0 


9 


ARS-BFGL-NGS-23583 




132049236 


A4GNT 


14386 


6.09E-06 


0.247 


CI 0:0 


15 


BTB-00059412 




132069357 


A4GNT 


1094 


1.62E-05 


0.658 


CI 0:0 


10 


ARS-BFGL-NGS-37095 




132497074 


S0X14 


46530 


1 .23E-05 


0.500 


C10:0 


17 


BTA-51403-no-rs 




132518716 


S0X14 


68172 


1.91E-05 


0.774 


CI 0:0 


3 


ARS-BFGL-NGS-91327 


2 


12533969 


ZNF804A 


593727 


3.70E-06 


0.150 


CI 4:0 


27 


ARS-BFGL-NGS-91 327 


2 


12533969 


ZNF804A 


593727 


2.19E-05 


0.889 


CI 4:0 


17 


ARS-BFGL-NGS-1 17409 


2 


37256390 


TANO 


Within 


7.04E-06 


0.286 


CI 4:0 


10 


Hapmap42557-BTA-47352 


2 


38748216 


LOCI 01 907729 


Within 


5.01 E-06 


0.204 


CI 4:0 


8 


BTB-01 053755 


2 


40313359 


NR4A2 


295672 


4.28E-06 


0.174 


C12:0 


21 


ARS-BFGL-NGS-30621 


2 


131539751 


HSPG2 


Within 


2.1 lE-05 


0.855 


CI 2:0 


15 


ARS-BFGL-NGS-20205 


3 


92427474 


SSBP3 


Within 


1 .40E-05 


0.569 


CI 4:0 


20 


BTB-01 477571 


5 


43481128 


CN0T2 


Within 


1 .25E-05 


0.509 


CI 4:0 


1 


Hapmap49848-BTA-1 06779 


5 


45089737 


CPM 


Within 


1.58E-07 


e.42E-03 


CI 0:0 


6 


Hapmap49071 -BTA-1 7699 


5 


92618397 


PIK3C2G 


1 15016 


7.53E-06 


0.306 


CI 2:0 


10 


BTB-01019973 


7 


79747454 


LOC101904982 


440947 


7.57E-06 


0.308 


C12:0 


7 


ARS-BFGL-NGS-67383 


7 


108307729 


EFNA5 


738306 


6.00E-06 


0.244 


CI 0:0 


1 


BTB-01 5561 97 


9 


16892513 


HTRIB 


409905 


5.89E-07 


2.39E-02 


CI 2:0 


3 


BTB-01 5561 97 


9 


16892513 


HTR1B 


409905 


2.66E-06 


0.108 


CI 2:0 


20 


ARS-BFGL-NGS-1 0471 9 


9 


18955509 


HMGN3 


44851 


2.10E-05 


0.851 


C10:0 


2 


ARS-BFGL-BAC-35400 


9 


21 165167 


FAM46A 


594573 


2.28E-06 


0.093 


CI 4:0 


18 


ARS-BFGL-BAC-35400 


9 


21 165167 


FAM46A 


594573 


7.23E-06 


0.294 


C10:0 


5 


ARS-BFGL-NGS-61979 


9 


23001645 


UBE3D 


42877 


6.06E-06 


0.246 


CI 2:0 


6 


ARS-BFGL-NGS-61979 


9 


23001645 


UBE3D 


42877 


5.66E-06 


0.230 


CI 0:0 


9 


Hapmap39984-BTA-21408 


9 


28538817 


LOCI 00848869 


46056 


1.10E-05 


0.447 


CI 4:0 


25 


Hapmap41 1 09-BTA-93077 


11 


42713681 


BCLUA 


358295 


1 .83E-05 


0.743 


CI 4:0 


5 


ARS-BFGL-BAC-5848 


12 


68657690 


GPC6 


Within 


2.29E-06 


0.093 


CI 0:0 


20 


BTA-37592-no-rs 


15 


72716471 


LRRC4C 


35911 


2.17E-05 


0.881 


CI 4:0 


23 


BTB-00634528 


16 


31813761 


SMYD3 


Within 


1.35E-05 


0.549 


CI 4:0 


6 


UA-IFASA-8132 


16 


33607353 


C16Hlorfl00 


7625 


2.73E-06 


0.111 


C10:0 


21 


BTB-00648332 


16 


55421856 


PDPN 


Within 


2.40E-05 


0.973 


CI 4:0 


16 


ARS-BFGL-NGS-1 02640 


17 


5963196 


PETl 12 


17972 


6.63 E-06 


0.269 


CI 4:0 


2 


BTB-00669395 


17 


6266432 


FAM160A1 


Within 


2.09E-06 


0.085 


CI 4:0 


4 


Hapmap47945-BTA-41 852 


17 


6295259 


FAM160A1 


Within 


2.26E-06 


0.092 


CI 4:0 


3 


BTB-00669586 


17 


6322271 


FAM160A1 


Within 


2.09E-06 


0.085 


CI 4:0 


13 


ARS-BFGL-NGS-20893 


17 


6669905 


SH3D19 


Within 


6.41 E-06 


0.260 


CI 4:0 


14 


ARS-BFGL-NGS-1 15234 


18 


13136171 


JPH3 


49852 


6.48E-06 


0.263 


CI 0:0 


14 


Hapmap23685-BTA-1 32541 


18 


54271729 


STRN4/SPHK2 


Within/1442497 


1 .59E-05 


0.647 


C10:0 


12 


UA-IFASA-7471 


18 


54311149 


SLC1A5/SPHK2 


Within/1403077 


1.41E-05 


0.571 


CI 0:0 


19 


ARS-BFGL-NGS-34500 


19 


15710458 


TMEM132E/ACACA 


81903/1714400 


2.11E-05 


0.857 


CI 0:0 


7 


ARS-BFGL-NGS-39328 


19 


51326750 


CCDC57/FASN 


Within/58172 


8.54E-06 


0.347 


CI 2:0 


2 


ARS-BFGL-NGS-39328 


19 


51326750 


CCDC57/FASN 


Within/58172 


l.ieE-06 


4.71 E-02 


CI 4:0 


11 


ARS-BFGL-NGS-39328 


19 


51326750 


CCDC57/FASN 


Within/58172 


6.01 E-06 


0.244 


CI 4:0 


21 


ARS-BFGL-NGS-87102 


20 


49859323 


CDH12 


585652 


1 .26E-05 


0.513 


CI 4:0 


12 


ARS-BFGL-NGS-1 11676 


20 


51073910 


CDH12 


Within 


6.02E-06 


0.244 


CI 2:0 


13 


Hapmap53927-rs29025287 


20 


53303717 


CDH18 


106159 


1.03E-05 


0.417 



PLOS ONE I www.plosone.org 3 May 2014 | Volume 9 | Issue 5 | e96186 



GWAS for Milk Fatty Acid Traits in Holstein 



Table 2. Cont. 







Trait 


Rank° 


SNP name 


Chr. 


Positlon(bp) 


Nearest Gene/ 
Candidate Gene 


Distance 
(bp) 


Raw P_value 


P_value 
Bonferroni 


C12:0 


11 


BTB-00787949 


20 


53333822 


CDH18 


76054 


8.44E-06 


0.343 


CI 4:0 


9 


BTB-01 583562 


20 


55425112 


LOC784462 


71405 


4.71E-06 


0.191 


C10:0 


16 


BTA-12468-no-rs 


21 


9375095 


ARRDC4/IGF1R 


436067/1107002 


1 .85E-05 


0.749 


CI 2:0 


12 


BTA-12468-no-rs 


21 


9375095 


ARRDC4/ICF1R 


436067/1107002 


l.OlE-05 


0.411 


CI 0:0 


8 


BTA-76414-no-rs 


21 


9528223 


ARRDC4/IGF1R 


589195/1260130 


9.90E-06 


0.402 


CI 2:0 


1 


BTA-76414-no-rs 


21 


9528223 


ARRDC4/IGF1R 


589195/1260130 


3.94E-07 


1.60E-02 


CI 0:0 


18 


ARS-BFGL-NGS-40159 


21 


21142616 


FANCVRLBPl/PLINl 


Within/Within/ 
360208 


1 .94E-05 


0.786 


CI 2:0 


14 


Hapmap26394-BTA-1 36497 


22 


27309084 


CNTN3 


Within 


1 .05E-05 


0.424 


CI 4:0 


7 


napmap2D394-D 1 A-1 3d49/ 


22 


z/3UyUo4 


CNTN3 


Within 


3.50E-06 


0.142 


CI 0:0 


1 1 


HapmapD /060-rs29023D 1 0 


24 


4805/59 


FBX015 


246097 


1 .33E-05 


0.539 


CI 2:0 


1 9 


napmapD/06U-rs29023Dl 0 


24 


^ O AC "7 C A 

4805/59 


FBX015 


246097 


1 .96E-05 


0.795 


CI 2:0 


4 


ADC DC/~I M/~C ~JOAr\~I 

AKb-DrLiL-NLib-/849/ 


24 


21330516 


SLC39A6 


Within 


T OAC AC 

2.8yb-06 


0.1 1 7 


CI 2:0 


8 


DTD nnOOCCIT 

D 1 D-UUbeSDD 1 Z 


24 


^AATQ-y-yc 
3UU2e5//5 


CHST9 


831 4 


6.08 E-06 


0.247 


CI 0:0 


13 


BTB-01 077939 


26 


7685110 


PRKGl 


Within 


1 .44E-05 


0.586 


C12:0 


22 


BTA-111275-no-rs 


26 


9195089 


MINPPl/PRKGl 


Within/851454 


2.39E-05 


0.969 


CI 4:0 


22 


ARS-BFGL-NGS-1 13226 


27 


97306 


LOa 00335608 


Within 


1 .34E-05 


0.544 


CI 2:0 


15 


BTB-01 926888 


27 


16398882 


TRML2/ACSL1 


303782/2110549 


1 .34E-05 


0.545 


CI 2:0 


5 


BTB-01 603522 


27 


16421445 


JRML2/ACSL1 


281219/2133112 


3.55E-06 


0.144 


C12:0 


18 


ARS-BFGL-NGS-1 8922 


29 


21930571 


LUZP2 


1125590 


1 .72E-05 


0.696 


CI 4:0 


24 


ARS-BFGL-NGS-1 9057 


29 


44196154 


CDC42EP2 


1060 


1.58E-05 


0.640 


CI 4:0 


15 


Hapmap60349-r529021 239 


X 


14062133 


ZNF280C 


Within 


6.51 E-06 


0.264 



Note: 

*^Ranl< represents ranking of significant SNPs within each oftrait; The P_value with bold type represents the significance of genome-wise level; The gene name with bold 
type represents the nearest known gene to the significant SNPs; The gene name with bold type and underline represents the nearest novel candidate gene to the 
significant SNPs. 

doi:1 0.1 371 /journal.pone.00961 86.t002 



mutations is still a challenge due to the commonly existing 
limitations [15]. 

At present, genome-wide association study (GWAS) has become 
a powerful strategy to identify genetic variants associated with 
complex traits. Since the first GWAS was published in 2005 [16], a 
great number of relative studies were conducted in human and 
domestic animals. Of them, several GWASs have been applied to 
detect genes or loci for milk production traits [17,18,19], 
conformation traits [19], reproduction traits [20,21], healthy traits 
[22,23], etc, in dairy and beef cattle. However, only studies have 
been carried out for fatty acids in Dutch dairy cattle [24,25] . We 
herein performed a GWAS for 22 milk fatty acid traits in a 
Chinese Holstein population to identify genes and chromosome 
segments with large effects on such traits. 

Materials and Methods 

The milk samples were collected along the regxilar quarantine 
inspection of the farms. The whole procedure for sample collection 
was carried out in strict accordance with the protocol approved by 
the Animal Welfare Committee of China Agricultural University 
(Permit Number: DK996). 

Phenotypic data and traits 

The Chinese Holstein population in this study comprised 784 
cows, the daughters of 21 sire families. AH cows in this study were 



from 18 farms of the Beijing Sanyuan Dairy Farm Center, where 
routine standard performance test, i.e. Dairy Herd Improvement 
system (DHI) have been carried out since 1999. A total of 50 ml 
milk sample was collected for each cow from the DHI laboratory 
of the Beijing Dairy Cattle Center during November to December, 
2012. The procedure of milk sample collection was carried out 
corresponding to DHI sampling (dairy herd improvement). After 
DHI measure, the remaining milk samples were taken back to the 
laboratory within 4°C cooler and then stored at — 20°C. 

Phenotypic values of 16 kinds of main milk fatty acids were 
measured by gas chromatography at the Ministry of Agriculture 
Feed Industry Centre of China (http://www.mafic.ac.cn/intro/ 
default.asp), which included SFA of C8:0, C10:0, C12:0, C14:0, 
C16:0, C18:0, C20:0, C22:0; MUFA of C14:l, C16:l, C18:ln9c; 
PUFA of CLA (cM-9, trans-ll C18:2), C18:3n3, C18:3n6, 
C18:2n6c and C20:5n3. Before measuring, milk samples should 
be done with pretreatment. First, total milk fat were extracted 
from approximately 2 ml of each milk sample. The specific 
procedure was as follows: 2 ml milk was mixed with 4 ml solution 
of N-hexane/isopropyl alcohol (3:2) and 2 ml solution of NajSO^, 
and centrifuged at 3,000 for 20 min. The upper layer was 
collected into 20 ml hydrolysis tube and 200 pi of C19:0 methyl 
ester as the internal standard was mixed, and then the extracted fat 
was dried under nitrogen. Methyl esters of fat were performed in 
the next step. 2 ml of NaOCH3/Methanol was put into the above 
hydrolysis tube for 15 min water bath under 50°C, and was mixed 



PLOS ONE I www.plosone.org 



4 



May 2014 | Volume 9 | Issue 5 | e96186 



GWAS for Milk Fatty Acid Traits in Holstein 
Table 3. Genome-wise and suggestive significant SNPs for long-chain saturated fatty acid traits (LCFA). 



Nearest Gene/ Raw P_value 



Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Candidate Gene 


Distance(bp) 


P_value" 


Bonferroni 


CI 6:0 


1 


ARS-BFGL-NGS-1 11692 


0 


0 


NA 


NA 


1.84E-06 


0.075 


C22:0 


1 


ARS-BFGL-NGS-1 09692 


1 


15796320 


NCAM2 


432009 


6.70E-07 


2.72E-02 


CI 8:0 


78 


ARS-BFGL-NGS-1 11491 


1 


68457522 


PTPLB 


6964 


1 .46E-05 


0.593 


CI 8:0 


91 


Hapmap50666-BTA-34589 


1 


68533156 


PTPLB 


Within 


1 .88E-05 


0.761 


CI 8:0 


8 


ARS-BFGL-NGS-761 1 1 


1 


103219606 


SI 


14878 


7.51 E-07 


3.05E-02 


CI 8:0 


25 


BTB-00048807 


1 


106245603 


OTOU 


454829 


2.46E-06 


0.100 


CI 8:0 


26 


BTB-00048739 


1 


106295981 


OTOLl 


404451 


2.77E-06 


0.112 


CI 8:0 


16 


ARS-BFGL-NGS-111111 


1 


146302724 


HSF2BP/AGPAT3 


43214/402889 


1.68E-06 


0.068 


CI 8:0 


41 


ARS-BFGL-NGS-1 09493 


1 


146354654 


HSF2BP/ACPAT3 


Within/350959 


5.28E-06 


0.214 


C18:0 


42 


BTA-56389-no-rs 


1 


146384457 


HSF2BP/AGPAT3 


Within/321156 


5.28E-06 


0.214 


CI 8:0 


93 


Haprinap5991 7-rs2901 241 8 


2 


24519348 


METAPID 


Within 


1.92E-05 


0.778 


CI 8:0 


88 


ARS-BFGL-BAC-2793 


2 


47449015 


KIF5C 


Within 


1 .79E-05 


0.726 


CI 8:0 


66 


Hapmap53388-rs2901 0903 


2 


63581955 


MGAT5 


129790 


1.18E-05 


0.479 


C18:0 


63 


BTB-01373917 


2 


79056259 


GYPC/STATl 


199326/837973 


1.12E-05 


0.453 


CI 8'0 


22 




2 




\3 r rL/ J 1 M 1 1 


ojODd/ DUO 1 'ty 


2 26E-06 


0 092 


CI 8:0 


52 


ADC DP("I M("C QQn^n 
Anj>-DroL-INU3-yyU3U 


2 


QQ 1 en 1 Q1 

yo 1 DU 1 y 1 


1 Shiran 


\A/i-i-l-t 1 1-1 

vviinin 


y. 1 4t-uo 


0,371 


CI 8:0 


1 1 


ADC RPf~l Mf-Q /IC^^QI 
Anb-bruL-INuD-^Doy 1 


2 


1 zo4e54/yu 


D{ IMV7/CADDJ 


i^^tSoO/D/uuyou 


9.56E-07 


3.88E-02 


CI 8'0 


30 




2 


1 zoDzy 1 uz 


nUIMAj/ rMDrJ 


\ UUD / 'f/ D / 4 J Z / Z 


3 22E-06 


0 131 


CI 8'0 


1 Q4 


ADC RPf^l Mr^C CROC; 


2 


1 jjOZU I // 


HTR6 


vvitnin 


2 41 E-05 


0 978 


CI 8'0 


1 0 


ADC Rpr;! Mr^c ^icnn^ 


2 




IFF02 


46005 


8 99E-07 


3 65E-02 


CI 8'0 


50 


ADC RFf^l Mf:;C 1 CflRT 




991 777 


MPZL 1 


vviimn 


H 7nF C\fi 
O. / ut uo 


0 353 


C18'0 


55 


napmapj^is J J 
BES3_Contig373_1182 




5549277 


IKZFl 


31 842 


1 OOE-05 


0 407 


C18:0 


79 


BTA-122414-no-rs 


4 


34967013 


SEMA3D 


596525 


1 .50E-05 


0.610 


C16:0 


12 


Hapmap40292-BTA-71 565 


4 


81400732 


C4H7orflO 


Within 


1 .39E-05 


0.564 


CI 8:0 


7 


Hapmap30257-BTA-1 42970 


5 


25358659 


USP44 


44459 


6.91 E-07 


2.81 E-02 


CI 8:0 


60 


Hapmap41951-BTA-73168 


5 


28442563 


SLC4A8 


25682 


1 .05E-05 


0.425 


CI 8:0 


73 


Hapmap49848-BTA-1 06779 


5 


45089737 


CPM 


Within 


1.31 E-05 


0.533 


CI 8:0 


12 


Hapmap50366-BTA-46960 


5 


68610818 


CHSTl 1 


Within 


1 .02E-06 


4.14E-02 


CI 6:0 


3 


Hapmap39862-BTA-74478 


5 


85672503 


BCATl 


74712 


4.88E-06 


0.198 


C16:0 


7 


BTB-003 16348 


7 


64939808 


ATOXl 


Within 


1.21 E-05 


0.490 


CI 8:0 


82 


BTB-01 553821 


7 


107940320 


EFNA5 


1105715 


1.63E-05 


0.663 


C18:0 


75 


BTB-00995040 


8 


22411464 


MIR31 


123333 


1.32E-05 


0.536 


CI 8:0 


18 


BTB-01 709624 


8 


33060034 


LOC101904752 


1103805 


1 .86E-06 


0.075 


CI 8:0 


43 


BTA-99986-no-rs 


8 


33175156 


LOC101904752 


988683 


5.54E-06 


0.225 


CI 8:0 


61 


BTB-01 973796 


8 


34146856 


LOC101904752 


16983 


1.07E-05 


0.434 


CI 8:0 


24 


BTB-01 929442 


8 


34560427 


LOC101904752 


52452 


2.45E-06 


0.099 


CI 8:0 


70 


ARS-BFGL-NGS-1 7346 


8 


70653923 


PEBP4 


Within 


1.27E-05 


0.514 


CI 8:0 


67 


Hapmap57757-ss4652621 5 


8 


74949334 


BNIP3L 


Within 


1.18E-05 


0.481 


CI 8:0 


46 


Anj-t)ruL-|\u J- 1 KJD/ DO 


8 


/oucyj/ 1 


&DTV 
Art A 


71 56 


6.77 E-06 


0.275 


CI 8'0 


53 


D 1 D UUzZ/ JO 1 


o 
0 


/oil yuuz 


DM^ I& 1 
Ui\r\Jr\ 1 


1 365 


9 23E-06 


0 375 


CI 8'0 


29 


RTR nn^'^Qi 1 1 

D 1 [3 UUjDy 1 1 Z 


g 


7fiiin^A'5'5 
/ j'fZZ 


NFXl 


vviinin 


3 1 4E-06 


0 1 28 


CI 8:0 


28 


ARS-BFGL-NGS-9052 


8 


78009328 


FRMD3 


181107 


3.10E-06 


0.126 


CI 8:0 


94 


BTB-01 9003 16 


8 


78226921 


UBQLNl 


Within 


1.93E-05 


0.783 


C18:0 


15 


Hapmap23947-BTA-1 5301 3 


8 


78347212 


GKAP1 


Within 


1 .66E-06 


0.067 


CI 8:0 


37 


Hapmap54400-rs29020952 


8 


79534851 


NTRK2 


Within 


4.82E-06 


0.196 


CI 8:0 


76 


Hapmap34874- 
BES3_Contig415_1312 


8 


79559282 


NTRK2 


Within 


1 .34E-05 


0.545 


CI 8:0 


81 


ARS-BFGL-NGS-1 01 844 


9 


1 3292880 


SLCUAS 


12122 


1 .59E-05 


0.644 


CI 8:0 


96 


ARS-BFGL-NGS-1 01 978 


9 


18123535 


HTRIB 


818548 


2.04E-05 


0.830 
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Table 3. Cont. 



Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value" 


P_value 
Bonferroni 


CI 8:0 


74 


ARS-BFGL-NGS-1 17379 


9 


97397264 


S0D2 


1895 


1.32E-05 


0.535 


CI 8:0 


45 


ARS-BFGL-NGS-66090 


10 


1 1 1 04820 


SERINC5 


Within 


5.91 E-06 


0.240 


CI 8:0 


39 


BTA-23031-no-rs 


10 


27606965 


OR4K13 


6376 


4.93E-06 


0.200 


C16:0 


5 


ARS-BFGL-NGS-23492 


10 


34637715 


RASGRPl 


356723 


1.13E-05 


0.460 


C16:0 


10 


ARS-BFGL-NGS-503 


10 


34705386 


RASGRPl 


424394 


1.35E-05 


0.547 


C16:0 


11 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC3SA6 


35850 


1.37E-05 


0.555 


CI 6:0 


17 


ARS-BFGL-NGS-1 4667 


10 


81481148 


GALNTLl 


Within 


1 .79E-05 


0.728 


CI 6:0 


8 


BTA-77299-no-rs 


10 


81549606 


SLC39A9 


Within 


1 .24E-05 


0.503 


CI 6:0 


16 


Hapmap44279-BTA-75297 


10 


84178157 


RGS6 


Within 


1.75E-05 


0.712 


C18:0 


102 


BTB-00471219 


11 


32709462 


NRXNl 


Within 


2.25E-05 


0.915 


CI 8:0 


103 


Hapmap25798-BTA-1 26388 


11 


32731961 


NRXNl 


Within 


2.40E-05 


0.976 


CI 8:0 


9 


HapmapSl 531 -BTA-98947 


11 


57522675 


REG3A 


880289 


8.78E-07 


3.S7E-02 


CI 8:0 


6 


BTB-01 328920 


11 


57639084 


LRRTM4 


792049 


S.79E-07 


2.35E-02 


C18:0 


40 


BTA-33625-no-rs 


11 


58306286 


LRRTM4 


124847 


5. 07 E-06 


0.206 


CI 8:0 


27 


Hapmap47549-BTA-25561 


11 


69506822 


LCLATl 


108007 


2.85E-06 


0.116 


CI 8:0 


97 


BTB-00866714 


12 


65771899 


LOC101907906 


250357 


2.09E-05 


0.850 


CI 8:0 


62 


ARS-BFGL-BAC-5848 


12 


68657690 


GPC6 


Within 


1.10E-05 


0.446 


CI 8:0 


100 


HapmapSl 198-BTA-27036 


12 


68703184 


GPC6 


Within 


2.15E-05 


0.873 


CI 8:0 


1 


ARS-BFGL-BAC-13788 


12 


69512332 


OCT 


Within 


9.17E-09 


3.72E-04 


C18:0 


65 


Hapmap2351 1 -BTA-1 1 9303 


12 


69768137 


S0X21 


24828 


1.14E-05 


0.463 


CI 8:0 


2 


ADC DCf~\ Klf~C ACIDCS 


1 2 


/02d04d/ 


ABCC4 


51 353 


1 .05E-08 


4.24C-04 


CI 8:0 


101 


ADC DC/~\ Klf~C Ar\"i~IO 

AKb-bSr(jL-N(j!5-4Uz/o 


1 3 


6/1 22063 


DLLAr/ACD52 


Wlthin/2280D43 


2.25E-05 


0.914 


C16:0 


4 


ADC DC/~I M/~c Ar\r>n 


14 


1801 1 16 


DGATl 


Within 


9.U4h-06 


0.367 


CI 8:0 


57 


DTD r\r\n c ~7 ^ 

D 1 d-UUdd/dJ2 


14 


24d4326d 


XKR4 


28554 


1 .03E-05 


0.416 


CI 8:0 


64 


ADC DCf~\ Klf~C 1 1 nmi 


14 


40616090 


LRIbrLUl 


144633 


1.12E-05 


0.456 


CI 8:0 


36 


ARS-BFGL-NGb-2025 


14 


65452336 


ZNF706 


1 03905 


4.72E-06 


0.1 92 


CI 8:0 


19 


Hapmap35102- 
BES3_Contig324_378 


14 


66523007 


RGS22 


Within 


1 .87 E-06 


0.076 


CI 8:0 


17 


ARS-BFGL-BAC-1991 


14 


82648300 


SNX 16 


491346 


1.73E-06 


0.070 


CI 8:0 


34 


ARS-BFGL-NGS-1 2066 


15 


27291268 


BUD13 


542743 


3.98E-06 


0.162 


CI 8:0 


54 


BTB-01 465034 


15 


49987359 


LOC784376 


1153 


9.65E-06 


0.392 


C16:0 


2 


ARS-BFGL-NGS-97658 


15 


68069900 


C15Hllorf74 


158748 


2.72E-06 


0.110 


C18:0 


92 


BTA-37923-no-rs 


15 


85011763 


GLB1L2 


Within 


1 .88E-05 


0.761 


CI 8:0 


3 


ARS-BFGL-NGS-68533 


16 


60615012 


LOG 101 902340 


35998 


2.71 E-07 


1.10E-02 


CI 8:0 


31 


ARS-BFGL-NGS-1 02798 


16 


61477887 


RALGPS2 


82055 


3.31 E-06 


0.134 


CI 8:0 


95 


BTB-00653808 


16 


62179726 


AXDNDl 


Within 


2.01 E-05 


0.817 


CI 8:0 


87 


ARS-BFGL-NGS-1 17800 


16 


62871926 


LHX4 


Within 


1 .75E-05 


0.711 


CI 8:0 


83 


ARS-BFGL-NGS-1 02835 


17 


3855180 


SFRP2 


17042 


1.67E-05 


0.678 


CI 8:0 


69 


ARS-BFGL-NGS-1 00229 


17 


4153159 


KIAA0922 


Within 


1.21 E-05 


0.493 


CI 8:0 


14 


ARS-BFGL-NGS-77485 


17 


4926550 


TMEM154 


39267 


1.47 E-06 


0.059 


CI 8:0 


38 


ARS-BFGL-NGS-1 18611 


17 


5357573 


F8XW7 


23715 


4.91 E-06 


0.199 


CI 8:0 


89 


ARS-BFGL-NGS-38059 


17 


5829384 


PETl 12 


151784 


1 .79E-05 


0.728 


CI 8:0 


84 


ARS-BFGL-NGS-1 14953 


17 


9753430 


NR3C2 


Within 


1 .70E-05 


0.688 


CI 8:0 


98 


ARS-BFGL-NGS-1 11 098 


17 


67193210 


PIWIL3 


27876 


2.10E-05 


0.851 


C18:0 


13 


ARS-BFGL-NGS-71116 


17 


68002540 


MY0188 


Within 


1 .09E-06 


4.41 E-02 


CI 8:0 


86 


ARS-BFGL-NGS-75816 


19 


36674728 


ANKRD40/SREBF7 


Within/1424056 


1 .70E-05 


0.692 


CI 8:0 


71 


BTB-01 790846 


20 


13263157 


SREKl 


128568 


1.27E-05 


0.516 


CI 8:0 


99 


Haprinap56230-r529025779 


20 


25535571 


NDUFS4 


12189 


2.14E-05 


0.869 


CI 8:0 


105 


ARS-BFGL-NGS-1 18998 


20 


32030332 


GHR/OXCTl 


Within/653664 


2.43E-05 


0.988 
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Table 3. Cont. 





Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value" 


P value 
Bonferroni 


CI 8:0 


48 


ARS-BFGL-NGS-87102 


20 


49859323 


CDH)2 


585652 


7.71 E-06 


0.313 


CI 8:0 


23 


BTB-01 475042 


20 


51845309 


CDH12 


220958 


2.29E-06 


0.093 


CI 8:0 


20 


ARS-BFGL-NGS-51112 


20 


52816141 


corns 


593735 


1.89E-06 


0.077 


C16:0 


18 


BTB-00818669 


21 


40502782 


PRKDl 


25243 


1 .95E-05 


0.790 


C16:0 


13 


BTB-01 240884 


21 


41484660 


G2E3 


128084 


1 .48E-05 


0.600 


C16:0 


9 


Hapmap24313-BTA-29957 


21 


51294895 


LRFN5 


476940 


1 .30E-05 


0.526 


CI 8:0 


58 


AKb-DFGL-DAC-467u7 


23 


5756655 


TINAG 


196853 


1 .03E-05 


0.419 


CI 6:0 


15 


BTA-56520-no-rs 


23 


32581285 


LOC53701 7/B TN1A1 


59797/1214677 


1.73E-05 


0.700 


CI 8:0 


33 


Hapmap58547-rs29023020 


23 


3591 1017 


LOCI 0084795 1/PRL 


Within/797267 


3.81 E-06 


0.155 


CI 8:0 


49 


BTA-58814-no-rs 


24 


5828401 


CBLN2 


207345 


8.20E-06 


0.333 


CI 8:0 


80 


BTB-00878928 


24 


5883369 


CBLN2 


26231 3 


1.51 E-05 


0.612 


CI 8:0 


59 


Hapmap34424- 
BES10_Contig566_926 


24 


6408329 


CBLN2 


787273 


1 .04E-05 


0.421 


C18:0 


4 


ARS-BFGL-NGS-1 09955 


24 


11625175 


CDH7 


411164 


3.29E-07 


1 .33E-02 


CI 8:0 


51 


BTA-24495-no-rs 


24 


1 5527942 


PIK3C3 


1 1 80676 


8.82E-06 


0.358 


CI 8:0 


56 


ARS-BFGL-NGS-13314 


24 


45775076 


SLC14A2 


Withlin 


1.01 E-05 


0.408 


CI 8:0 


35 


BTA-59652-no-rs 


25 


19315456 


LOCS24391 


Witilin 


4.37E-06 


0.178 


CI 6:0 


6 


BTB-01 61 91 01 


26 


16035322 


CYP2C87 


Withiin 


1.16E-05 


0.472 


CI 8:0 


72 


ARS-BFGL-NGS-12828 


26 


37018466 


GFRAl 


Witilin 


1 .29E-05 


0.524 


CI 8:0 


21 


ARS-BFGL-NGS-1 00468 


26 


38060272 


PDZD8 


39152 


2.10E-06 


0.085 


CI 8:0 


47 


ARS-BFGL-NGS-1 1 1 901 


26 


41183634 


WDR1 1 


28744 


7.62E-06 


0.309 


CI 8:0 


68 


UA-IFASA-5698 


26 


42673967 


HTRA1 


Witilin 


1.21 E-05 


0.491 


CI 8:0 


77 


Hapmap6081 0-rs2901 2623 


27 


1235541 1 


LOC100848735 


1 30088 


1 .39E-05 


0.565 


C20:0 


2 


BTB-00965 1 97 


27 


2681 3022 


NRGl 


Within 


1 .93E-05 


0.783 


C20"0 


] 


ARS-RFGI -NGS-77nn'? 


27 


3\JOOD\JO 1 


UNC5D 


Within 

V V 1 LI III 1 


3 76E-06 


0 1 53 


CI 8:0 


85 


Hapnnap35611- 
SCAFFOLDl 20249_1 7244 


27 


44005137 


ZNF385D 


5580 


1 .70E-05 


0.691 


C16:0 


14 


Hapmap2741 8-BTA-l 47969 


29 


23469370 


LOCS40991 


19676 


1 .66E-05 


0.675 


CI 8:0 


32 


ARS-BFGL-NGS-1 9057 


29 


44196154 


CDC42EP2 


1060 


3. 65 E-06 


0.148 


CI 8:0 


5 


ARS-BFGL-NGS-1 1681 


X 


3622615 


Gene desert 


NA 


3.74E-07 


1 .52E-02 


CI 8:0 


44 


Hapnnap48540-BTA-97806 


X 


8700029 


ODZ; 


Witilin 


5. 77 E-06 


0.234 


C18:0 


90 


BTB-01 492502 


X 


8724630 


ODZ; 


Witilin 


1 .84E-05 


0.746 



Note: see note to Tabie 2. 

doi:1 0.1 371 /journal.pone.00961 86.t003 



with 2 ml of hydrochloric acid/methanol solution (1:10) for 1.5 h 
water bath under 80°C. After the temperature fell to room 
temperature level, 3 ml of water and 6 ml of n-hexane were put 
into above hydrolysis tube, mixed, vortexed, and stratified. The 
upper layer was collected and dried under nitrogen, and finally 
dissolved in 1 ml of n-hexane. 1 ml methyl esters of fatty acids 
were prepared to be determined by gas chromatography using a 
gas chromatograph (6890N, Agilent) equipped with a flame- 
ionization detector and a high polar fused silica capillary column 
(SP™-2560, 100 mx0.25mm ID, 0.20 ixm film). About 1 |il 
sample was injected under the following gas chromatography 
conditions: Hehum was used as the carrier gas at a flow of 45 ml/ 
min. The split ratio was 100:1. The oven temperature was 
programmed at 100°C and held for 10 min, then increased to 
160°C at a rate of 6°C/min, held for 10 min, increased to 200°C 
with 5°C/min, held for 20 min, increased to 240°C at a rate of 
4°C/min and held for 12 min. Both the injector temperature and 
the detector temperature were set on 260°C. Individual fatty acids 



were identified and quantified by comparmg the methyl ester 

chromatograms of the milk fat samples with the chromatograms of 

pure fatty acids methyl ester standards (Supelco^'^ 37 Component 

FAME Mix), and were measured as the weight proportion of total 

fat weight (wt/wt%). 

Based on the phenotypes of 16 milk fatty acids, 6 additional 

traits were obtained including SFA, UFA, SFA/UFA (the ratio of 

SFA to UFA), C 1 4 index, C 1 6 index and C 1 8 index. The 3 indices 

cis — 9 unsaturated 

were calculated as —. 1 r x 100 [261. 

CIS — 9 unsaturated + saturated 

The descriptive statistics of these 22 fatty acid traits are 

presented in Table 1. Both SFA and UFA accounted for 

approximately 96% (wt/wt) of total fat. 

Genotypes and quality control 

The cows were genotyped using the lUumina BovineSNP50 
BeadChip (lUumina Inc., San Diego, CA, US), of which, some 
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Table 4. Genome-wise and suggestive significant SNPs for monounsaturated and polyunsaturated fatty acid traits (IVIUFA and 
PUFA). 





Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Dlstance(bp) 


Raw 

P_value° 


P_value 
Bonferroni 


C14:l 


34 


ARS-BFGL-NGS-15914 


0 


0 


NA 


NA 


1 .20E-06 


4.87E-02 


CLA 


2 


ARS-BFGL-NGS-43953 


1 


113855984 


GPR149 


288255 


1.31 E-05 


0.531 


C14:l 


39 


ARS-BFGL-BAC-42858 


2 


33401608 


LOCS34542 


231071 


2.02E-06 


0.082 


CLA 


4 


BTB-01 649999 


3 


48062099 


RWDD3 


336059 


1.65 E-05 


0.672 


C18:ln9c 


9 


Hapmap42304-BTA-73062 


5 


5049476 


KRRVOSBPL8 


14567/770655 


1.25 E-05 


0.509 


C18:1n9c 


19 


ARS-BFGL-NGS-8796 


5 


29095603 


LOCS10716 


Within 


2.25E-05 


0.915 


C18:1n9c 


8 


ARS-BFGL-NGS-69056 


5 


42285835 


CPNE8 


135184 


1.18E-05 


0.479 


C18:ln9c 


10 


Hapmap39862-BTA-74478 


5 


85672503 


BCATl 


74712 


1 .27E-05 


0.515 


C18:ln9c 


6 


ARS-BFGL-NGS-99256 


5 


104714350 


VWF/OLR7 


Within/4458707 


8.88E-06 


0.361 


C14:l 


64 


ARS-BFGL-NGS-1 10361 


7 


15782979 


COL5A3 


Within 


2.03E-05 


0.826 


C14:l 


61 


ARS-BFGL-NGS-1 04050 


7 


62839580 


CSNKlAl 


13832 


1 .48E-05 


0.600 


C14:l 


65 


Hapmap50476-BTA-79543 


7 


62889548 


CSNKlAl 


Within 


2.17E-05 


0.882 


C18:ln9c 


18 


BTB-003 16291 


7 


64892251 


SPARC 


Within 


2.15E-05 


0.874 


C18:ln9c 


12 


BTB-003 16348 


7 


64939808 


ATOXl 


5940 


1.32E-05 


0.534 


CLA 


6 


BTB-01 541 157 


7 


84582299 


ATP6AP1L 


15187 


1 .78E-05 


0.722 


C18:2n6c 


5 


BTA-19330-no-rs 


7 


93289032 


ARRDC3 


35938 


9.00E-06 


0.366 


CLA 


5 


Hapmap27874-BTA-146513 


9 


53871363 


GPR63 


53530 


1 .66E-05 


0.676 


C18:ln9c 


15 


Hapmap501 26-BTA-83733 


9 


55449737 


LOC101907134 


18763 


1 .49E-05 


0.606 


C16:l 


6 


ARS-BFGL-NGS-1 15094 


9 


93183596 


TIAM2 


Within 


1 .88E-05 


0.762 


C16:l 


4 


ARS-BFGL-NGS-57866 


9 


93246183 


TFS1M/TIAM2 


Within 


1 .34E-05 


0.545 


C18:ln9c 


7 


BTB-01 332998 


10 


73466092 


SLC38A6 


14392 


1 .02E-05 


0.414 


C18:ln9c 


5 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC38A6 


35850 


5.75E-06 


0.233 


C16:l 


5 


Hapmap381 87-BTA-l 05082 


12 


69216840 


OCT 


286101 


1.51 E-05 


0.614 


C16:l 


2 


Hapmap53988-rs29024591 


13 


65855988 


EPB41U/ACSS2 


23900/1014468 


9.31 E-06 


0.378 


C16:l 


3 


BPI-2 


13 


67833153 


L0CS14211/ACSS2 


Within/2991633 


1.16E-05 


0.470 


C18:2n6c 


4 


ARS-BFGL-NGS-57820 


14 


1651311 


LOC 10029491 6/ DCATI 


Within/144114 


2.47E-06 


0.100 


C14:l 


60 


ARS-BFGL-NGS-4939 


14 


1801116 


DCAT1 


Within 


1.45 E-05 


0.590 


C18:ln9c 


2 


ARS-BFGL-NGS-4939 


14 


1801116 


DGAT1 


Within 


3.01 E-06 


0.122 


C18:2n6c 


1 


ARS-BFGL-NGS-4939 


14 


1801116 


DGAT1 


Within 


6.87E-08 


2.79E-03 


C14:l 


25 


Hapmap52798-ss46526455 


14 


1923292 


MAF1/DGAT1 


Within/118454 


2.32E-07 


9.42E-03 


C18:2n6c 


3 


ARS-BFGL-NGS-1 07379 


14 


2054457 


LOC786966/DGAT1 


460/249619 


1.13E-06 


4.57E-02 


C18:ln9c 


1 


ARS-BFGL-NGS-1 00480 


14 


4364952 


TRAPPC9 


Within 


1.95 E-06 


0.079 


CI 4:1 


43 


Hapmap32234-BTC-0481 99 


14 


7314869 


KHDRBS3 


1 1 7233 


3.48E-06 


0.141 


C18:2n6c 


2 


ARS-BFGL-NGS-1 14448 


16 


32095335 


SMYD3 


Within 


4.72E-07 


1 .92E-02 


C16:l 


1 


BTB-01 090859 


16 


63248646 


XPRl 


51157 


4.29E-06 


0.174 


C14:l 


67 


ARS-BFGL-NGS-66923 


17 


18910914 


CCRN4L 


Within 


2.36E-05 


0.958 


C14:l 


37 


ARS-BFGL-NGS-1 02933 


17 


19315294 


CCRN4L 


404368 


1 .46E-06 


0.059 


CI 4:1 


16 


BTA-91575-no-rs 


17 


20099837 


SLC7A11 


291555 


3.42E-08 


1 .39E-03 


C14:l 


47 


ARS-BFGL-NGS-25840 


17 


23694684 


LOC783956 


163671 


4.48E-06 


0.182 


C14:l 


32 


BTB-01 585209 


17 


44910178 


ZNF605 


53460 


9.81 E-07 


3.98E-02 


CI 4:1 


17 


ARS-BFGL-NGS-1 09854 


17 


44963124 


ZNF605 


514 


4.23E-08 


1 .72E-03 


C18:ln9c 


16 


ARS-BFGL-NGS-87368 


19 


7762820 


C19H17orf67 


67914 


1 .76E-05 


0.713 


C14:l 


53 


BTA-117074-no-rs 


20 


247319 


PANK3 


42534 


7.72E-06 


0.313 


C14:l 


62 


ARS-BFGL-NGS-1 14602 


20 


1052840 


SLIT3 


979 


1.52E-05 


0.618 


C14:l 


57 


ARS-BFGL-NGS-1 01 925 


20 


1202954 


SLIT3 


151093 


8.96E-06 


0.364 


C18:ln9c 


4 


ARS-BFGL-NGS-1 16806 


20 


36450009 


GDNF 


180881 


4.66E-06 


0.189 


C18:ln9c 


13 


ARS-BFGL-NGS-14031 


20 


36561330 


GDNF 


69560 


1.46E-05 


0.594 



PLOS ONE I www.plosone.org 



8 



May 2014 | Volume 9 | Issue 5 | e96186 



GWAS for Milk Fatty Acid Traits in Holstein 



Table 4. Cont. 





Trait 


Rank^ 


SNP 


Chr. 


Posit ion (bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P value^ 


P_value 
Bonf erroni 


C18:ln9c 


1 1 


ARS-USMARC-Parent- 
DQ990835-rs29012811 


20 


36570529 


GDNF 


60361 


1 .28E-05 


0.519 


CLA 


1 


BTA-50525-no-rs 


20 


36917645 


WDR70 


Within 


1.16E-05 


0.469 


C18:ln9c 


20 


ARS-BFGL-NGS-1 11420 


20 


37553988 


SLC1A3 


42121 


2.30E-05 


0.935 


C18:ln9c 


17 


Hapmap57276-ss46526009 


20 


37596667 


SLC1A3 


Within 


2.04E-05 


0.829 


C18:2n6c 


6 


BTB-01 340958 


20 


60400902 


DNAH5 


840190 


1 .70E-05 


0.689 


CLA 


7 


Hapnnap2651 2-BTA-52638 


21 


55483469 


MIS18BP1 


175 


1.85E-05 


0.752 


C14:l 


54 


ARS-BFGL-BAC-46707 


23 


5756655 


TINAG 


196853 


7.82E-06 


0.318 


CLA 


3 


BTA-24495-no-rs 


24 


15527942 


LOC783699 


51814 


1.57E-05 


0.637 


C18:ln9c 


3 


ARS-BFGL-BAC-28144 


25 


2606575 


LOC788915 


21466 


4.34E-06 


0.176 


C14:l 


48 


Hapmap54064-rs2901 1 996 


26 


5526925 


PCDH15 


Within 


4.58E-06 


0.186 


C14:l 


35 


ARS-BFGL-NGS-1 3746 


26 


9866940 


RNLS 


65137 


1.42E-06 


0.058 


C14:l 


59 


BTA-61921-no-rs 


26 


10255258 


LIPJ 


7406 


1.42E-05 


0.577 


C14:l 


21 


ARS-BFGL-NGS-21794 


26 


10397362 


LOC100336557/LIPK 


16309/35288 


9.57E-08 


3.88E-03 


C14:l 


51 


BTA-08775-rs29022332 


26 


1 1201 198 


SLC16A12 


39720 


6.64E-06 


0.269 


C14:l 


44 


ARS-BFGL-NGS-53115 


26 


11528933 


KIF20B 


61222 


3.82E-06 


0.155 


C14:l 


24 


ARS-BFGL-NGS-63853 


26 


11942868 


MIR2895 


67662 


2.07E-07 


8.42E-03 


C14:l 


36 


ARS-BFGL-NGS-1 2381 


26 


12200948 


LOC100141242 


24477 


1.45E-06 


0.059 


C14:l 


20 


BTB-01 90841 7 


26 


12268427 


LOC784522 


15486 


5.62E-08 


2.28E-03 


C14:l 


66 


ARS-BFGL-NGS-41148 


26 


12364119 


HTR7 


10881 


2.20E-05 


0.892 


C14:l 


28 


Hapnnap5281 7-rs2901 1 969 


26 


14155229 


HHEX 


29160 


3.97E-07 


1.61E-02 


C14:l 


38 


ARS-BFGL-NGS-1 10475 


26 


15604631 


PLCEl 


Within 


1 .54E-06 


0.063 


CI 4:1 


19 


ARS-BFGL-NGS-29299 


26 


16614068 


PDUMl/SORBSl 


Within/73061 


5.42E-08 


2.20E-03 


C14:l 


23 


Hapmap41 595-BTA-60800 


26 


16791783 


SORBS! 


Within 


2.07E-07 


8.41 E-03 


C14:l 


40 


Hapmap58930-rs2901 0490 


26 


16822073 


SORBS 1 


Within 


2.34E-06 


0.095 


CI 4:1 


42 


ARS-BFGL-NGS-41056 


26 


18906121 


CRTACl 


Within 


2.70E-06 


0.110 


CI 4:1 


45 


ARS-BFGL-NGS-1 16902 


26 


18967997 


CRTACl 


Within 


4.31 E-06 


0.175 


C14:l 


55 


ARS-BFGL-NGS-25126 


26 


18994785 


CRTACl 


Within 


7.84E-06 


0.318 


C14:l 


15 


ARS-BFGL-NGS-23064 


26 


20365711 


NKX2-3/SCD7 


34516/767033 


2.93E-08 


1.19E-03 


C14:l 


14 


ARS-BFGL-NGS-77668 


26 


20393457 


NKX2-3/SCD7 


6770/744488 


2.23E-08 


9.05E-04 


C14:l 


33 


BTB-00930925 


26 


20474308 


SLC25A28/SCD1 


Within/658436 


1 .08E-06 


4.39E-02 


C14:l 


22 


ARS-BFGL-NGS-39397 


26 


20716721 


DNMBP/SCDl 


Within/416023 


1.37E-07 


5.55E-03 


C14:l 


27 


BTB-00930720 


26 


20903573 


LOC5U498/SCD1 


Within/244744 


2.8eE-07 


1.ieE-02 


C14:l 


56 


Hapmap4641 1-BTA-l 5820 


26 


20984335 


CHUK/SCDl 


Within/148409 


8.29E-06 


0.336 


C14:l 


31 


Hapmap31 825-BTA-l 58647 


26 


21056547 


PKD2U/SCD1 


Within/76197 


S.10E-07 


2.07E-02 


C14:l 


11 


Hapmap33073-BTA-1 62864 


26 


21180893 


SCD1 


32576 


1.31E-08 


5.33E-04 


C14:l 


1 


BTB-00931481 


26 


21226405 


WNT8B/SCD1 


14100/78088 


7.08E-13 


2.87E-08 


C14:l 


12 


ARS-BFGL-NGS-1 10077 


26 


21322557 


HIF1AN/SCDJ 


22399/174240 


1.54E-08 


e.23E-04 


CI 4:1 


10 


ARS-BFGL-NGS-1 08305 


26 


21363670 


HIFIAN/SCDJ 


63512/215353 


8.07E-09 


3.28E-04 


C14:l 


9 


BTB-a0931586 


26 


21409429 


PAX2/SCD7 


61334/261112 


2.25E-09 


9.12E-05 


C14:l 


5 


ARS-BFGL-NGS-1 14149 


26 


21702714 


LZTS2/SCD7 


656/564769 


5.50E-10 


2.23E-05 


CI 4:1 


6 


ARS-BFGL-NGS-1 16481 


26 


21977581 


LOC100847491/SCD1 


10062/829264 


5.77E-10 


2.34E-05 


C14:l 


7 


Hapnnap24832-BTA-1 38805 


26 


22016380 


BTRC/SCD1 


Within/868063 


7.18E-10 


2.92E-05 


C14:l 


8 


ARS-BFGL-NGS-6259 


26 


22059103 


BTRC/SCDl 


Within/910786 


7.23E-10 


2.93E-05 


CI 4:1 


4 


BTB-00932332 


26 


22118554 


BTRC/SCD1 


Within/970237 


3.56E-10 


1.45E-05 


C14:l 


3 


ARS-BFGL-NGS-1 07403 


26 


22889812 


NFKB2 


1586 


2.22E-12 


9.02E-08 


C14:l 


18 


Hapmap48222-BTA-1 22240 


26 


23641881 


C26HWorf26 


Within 


4.53E-08 


1.84E-03 


C14:l 


29 


Hapmap49372-BTA-91 009 


26 


23689229 


C26HWorf26 


3017 


4.e6E-07 


1 .89E-02 
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Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Dlstance(bp) 


Raw 

P.value" 


P_value 
Bonferroni 


C14:l 


46 


BTA-60918-no-rs 


26 


23853334 


CNNM2 


Within 


4.40E-06 


0.179 


C14:l 


13 


BTA-60935-no-rs 


26 


23876476 


CNNM2 


Witliin 


1.5eE-08 


e.33E-04 


C14:l 


52 


ARS-BFGL-NGS-2180 


26 


24477962 


SH3PXD2A 


Within 


7.46E-06 


0.303 


C14:l 


49 


ARS-BFGL-NGS-1092 


26 


24531763 


SH3PXD2A 


Within 


4.90E-06 


0.199 


CI 4:1 


2 


ARS-BFGL-NGS-118189 


26 


24786731 


SLK 


Within 


1.69E-12 


6.87E-08 


C14:l 


58 


UA-IFASA-4715 


26 


25314352 


CCDC147 


27285 


9.16E-06 


0.372 


C14:l 


30 


Hapmap28763-BTA-1 62328 


26 


26757136 


S0RCS3/ECHS1 


358845/891528 


5.05E-07 


2.05E-02 


C14:l 


26 


BTA-87355-no-rs 


26 


27251857 


SORCSl 


558341 


2.59E-07 


1 .05E-02 


C14:l 


41 


BTA-10873-rs29016424 


28 


18317414 


RTKN2 


41359 


2.35E-06 


0.095 


C14:l 


63 


ARS-BFGL-NGS-1 18293 


28 


45409471 


CXCLU 


Within 


1 .77E-05 


0.718 


C18:ln9c 


14 


BTB-01088519 


29 


22161137 


CCDC179 


85675 


1 .49E-05 


0.605 


C14:l 


50 


ARS-BFGL-NGS-27560 


X 


146596721 


STS 


Within 


5.19E-06 


0.211 



Note: see note to Table 2. 

doi:l 0.1 371 /journal.pone.00961 86.t004 



individuals were genotyped with die 54K cliip version 1 containing 
54,001 SNPs, and others were genotyped with the 54K chip 
version 2 including 54,609 SNPs. After genotype imputation by 
BEAGLE software (http://faculty.washington.edu/browning/ 
beagle/beagle.htnil), the common SNP markers in both version 
chips were used in this study, as a result, the total number of SNPs 
in the panel was 52,340. The SNP positions were based on the 
bovinegenome assembly UMD_3.1.66 (http://www.ncbi.nlm.nih. 
gov/genome/guide/ cow/). 

The quality control procedure was as follows, 20 daughters were 
excluded due to low call rate (<90%), leading to 764 daughters 
remaining for the association analysis. On the other hand, 1 1,736 
SNPs were removed for falling to meet the following requirements: 
652 SNPs with <90% genotype call rate, 10,798 SNPs with a 
minor allele frequency (MAF) <0.05, 286 SNPs with extreme 
value of Hardy-Weinberg equilibrium statistics (P<10 ). Even- 
tually, 40,604 SNPs passed these quality control filters, which was 
77.6% of the SNPs in the panel. The average distance between 
adjacent markers was quite constant among different chromo- 
somes. The shortest average distance was 56 kb on BTA25, and 
the longest average distance was 75 kb on BTA5 (except for 
198 kb on BTAX). 

Statistical analysis 

The statistical tests followed a two-step analysis. For the first 
step, phenotypic values were corrected for frxed non-genetic effects 
by using SAS 9.1 general linear model (GLM) procedure. The 
statistical model was: yijki = ti + Fi + Pj +Lk + eijki, where jiju was 
the unadjusted phenotype; fi was the overall mean; was the fixed 
effect of farm; Pj was the fixed effect of parity; was the fixed 
effect of stage of lactation; ejju was the random residual. In the 
second step, genome-wide association analyses were performed 
with quantitative trait procedure (additive model)of the PLINK 
software (vl.07) [27], and empirical p-values estimated based on 
the Wald-statistic. Individual pedigree of three generations was 
applied. 

Manhattan plots of genome-wide association analyses were 
produced with R2.15.1 software (http://www.r-project.org/). 



Significance level 

Bonferroni correction was apphed to adjust for multiple testing 
from the number of SNPs detected. A significant SNP at the 
genome-wise significance level was declared if a raw P value 
(unadjusted)<0.05/jV, jVis the number of SNP markers tested in 
analyses [28]. In the present study, Bonferroni genome-wise 
significance was 1.23E-06 (0.05/40604). As the Bonferroni 
correction threshold levels were strict and may lead to high false 
negatives, we calculated suggestive significant association threshold 
P-value as previously described [29], which was 2.46E-05 (1/ 
40604). 

Results 

The global view of P-values for all SNPs of each trait was shown 
in Additional file 1. In total, 83 genome-wise significant SNPs (P< 
1.23E-06) and 314 suggestive significant SNPs (P<2.46E-05) were 
detected for 22 milk fatty acids on all chromosomes, ranged from 3 
on BTA3, 4, 22 to 119 on BTA26 (Tables 2-6). For most of die 
studied fatty acids, significant associations were detected with more 
than one SNP, especially C14:l (67 SNPs: 34 genome-wide, 33 
suggestive), C18:0 (105 SNPs: 13 genome-wide, 92 suggestive), 
C14 index (84 SNPs: 49 genome-wide, 35 suggestive) and C18 
index (93 SNPs: 14 genome-wide, 79 suggestive) (Table 7). The 
most significant SNP (BTB-00931481) was associated with both 
C14 index (P=6.91E-17) and C14:l (P=7.08E-13) on BTA26. 
The top one common significant SNP (ARS-BFGL-NGS-4783) 
was associated with SFA, SFA/UFA and UFA. Besides, 146 SNPs 
were associated with multiple traits, especially ARS-BFGL-NGS- 
4939 on BTA14 for 9 traits. Further details on these associations 
are described as follows. 

Short- and medium-chain saturated fatty acid traits 
(SCFA and MCFA) 

For C10:0, C12:0 and C14:0, 21, 22 and 27 SNPs were 
detected, respectively. Of these 70 SNPs, 10 were associated with 
two or three traits. The most significant association of CI 0:0 
(P=5.89E-07), C12:0 (P= 3.94E-07), and C14:0 (P= 1.58E-07) 
were identified with BTB-01556197 on BTA9, BTA-76414-no-rs 
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Table 5. Genome-wise and suggestive significant SNPs for indices of fatty acid traits. 



Trait 


Rank^ 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value^ 


P_ value 
Bonferroni 


C 


4index 


26 


AR<;-RF(^l -NI(^<^-1 ^Q14 


Q 


Q 


NA 


NA 


2.59E-09 


1 .05E-04 


C 


4index 


38 


1 1 A IPAC A crot:') 
UA-lrAjA-DtSDZ 


0 


0 


MA 


MA 


O.U4C-UO 


3.26E-03 


c 


Sindex 


65 


D 1 D-UUZ/U 1 30 


0 


0 


MA 


MA 


1 .1 9E-05 


0.4S2 


c 


SIndex 


91 


l-l3nmar>ftnftA7-rc9Qn'57^41 
napi 1 idpouD^/ i iz:7UZ / jh- i 


Q 


Q 


NA 


NA 


2 25E-05 


0 912 


c 


Sindex 


23 


RTA-1 1 1771-nn-r<; 




37 T 74447 


EPHA3 


50365 


T RfiF-Ofi 
Z.OL;^ \j\J 


0 116 


c 


Sindex 


22 






41 577955 


EPHA6 


Within 

V V 1 LI III 1 


2 74E-06 


0 111 


c 


Sindex 


45 


/An J DrVJL INVJ J 1 1 J / Oj 




DO 1 Ji'D J J 


GTF2E1 


67771 


D. 1 UO 


0 249 


c 


Sindex 


35 


D 1 D-UUUJZZUU 




D/ / 044/0 




vviinin 


3.83E-06 


0.1 56 


Q 


Sindex 


33 






1 uzjo:?uu:? 


BCHE 


1 78966 


3 69E-06 


0 1 50 


Q 


Sindex 


24 


ARS-BFGL-NGS-111111 




146302724 


HSF2BP/AGPAT3 


43214/4028S9 


2.92E-06 


0.118 


Q 


Sindex 


12 


ARS-BFGL-NGS-1 09493 




146354654 


HSF2BP/AGPAT3 


Within/350959 


1.13E-06 


4.59E-02 


Q 


Sindex 


13 


BTA-56389-no-rs 




1463S4457 


HSF2BP/AGPAT3 


Within/321156 


1.13E-06 


4.59E-02 


Q 


Sindex 


32 


ARS-BFGL-NGS-76347 




146704618 


AGP AT 3 


995 


3.66E-06 


0.149 


Q 


Sindex 


84 


Hapmap5991 7-rs2901 241 8 


2 


24519348 


METAP1D 


Within 


1 .90E-05 


0.771 


Q 


SinoGX 


21 


Hapmap533S8-rs2901 0903 


2 


63581955 


MGAT5 


129790 


2.65 E-06 


0.107 


Q 


Sindex 


3 


ARS-BFGL-NGS-33744 


2 


7938S083 


GYPC/STAT1 


83687/506149 


8.95E-08 


3.63E-03 




Sinclcx 


31 


Hapmap5341 9-rs2901 51 59 


2 


88205436 


SAT 82 


43199 


3.57E-06 


0.145 




Sindsx 


5 


Hapmap33966- 
RFS? rnntln^ftR 774 


2 


88545567 


SATB2 


9001 S 


2.56E-07 


1 .04E-02 


C 


Sindex 


50 


AR<;-RF(^I -MC^-'>^R7'> 


2 


951 99905 


ADAM?? 


Within 

vv III III 1 


7 34E-06 


0 298 


c 


Sindex 


76 


ARc;-RF(^l -Mf^c;-QR^'^4 
/An J DrUL ino J iJO J jm- 


2 


9551 2347 


r/l J / AL/Z 


Within 

vv 1 11 III 1 


1 47E-05 


0 598 


c 


Sindex 


55 


AR<;-RF(^i -Mr;<^-QQn^n 


2 


1 DU 1 ^ 1 


/ INCRf) 


Within 
vv III III 1 


O.OJC UD 


0 359 


c 


Sindex 


1 5 


AR<;-RF(^I -M(^*^-imfiQ1 
/Anj DrVJL INU J 1 


2 


1 ZO^OH-/ ?U 


r\UI\AJ/ rrlDr J 


^AARRf\/^7^^Qf^^ 


141 E-06 


0 057 


c 


Sindex 


g 


AR<;-RF(^I 

MnJ DrUL INU J 1 1 O^ZH 


2 


1 ZODZ7 1 UZ 


Rl IMY^/FARP^ 

r\\Jt\A.J/ mOr J 


1 r\r\'^7A/'^7A^171 

1 UU J / H-/ J / ^ JZ / Z 


6.04E-07 


2.45E-02 


c 


Sindex 


28 


l-lanr-n;4t-i^n9^7-RTA-1 <19Q7n 
ndpi 1 IdpOUZD / D 1 r\ 1 H-Z3/ u 




zj J jod j:? 


U5P44 


44459 


3 24E-06 


0 131 


c 


Sindex 


75 


AR^^-RFni -Mr;<;-^Rn^R 

/An J DrVJL INVJ J jOU jO 




279921 79 


NR4A1 


Within 
vv III III 1 


1 45E-05 


0 5S8 


c 


Sindex 


59 


l-lanm3n419S1-RTA-731 68 


5 


ZOt-^tZ-JU J 


SLC4A8 


25682 


9 43 E-06 


0 3S3 


c 


Sindex 


16 


AR*;-RF(^I -MG*;-R7Qfi 


5 


99nQSfin3 

Z_7\JI?_JUUj 




Within 

V V 1 LI III 1 


1 82 E-06 


0 074 


c 


Sindex 


92 


AR^^-RFni -Mr;<^-'^^4RR 

/An J DrVJL InVJ J J jH-OO 




41 1 ^A^9R 


pRirwi F1 

r n/l_Ai_C / 


1 6691 5 


2 35E-05 


0 953 


c 


Sindex 


47 


AR*;-RF(^I -Mr;<^-1 1fiRQ7 
/Anj DrUL INU J 1 1 do:?/ 






PI RDI/ni R1 


WithinM'^nn'^7'^ 

vv ILI III 1/^ JUU JZ J 


ft AQF-flft 


0 264 


c 


Sindex 


73 


D 1 D U 1 DOJZ Ji? 




1 ZZOZOO 1 


UGT8 


81 40S 


1 44E-05 


0 5S6 


c 


Sindex 


1 4 


RTR-nn9dfii "^n 

□ 1 D VjUZH-D 1 D\J 




IHQQ^AIA 

ZU7:?3H-ZH 


PPA2 


Within 
vv ILI III 1 


1 .1 5E-06 


4.68E-02 


c 


Sindex 


30 


ndpi 1 idpzDUU 1 D 1 u joo i j 




AAQ7fi9A^ 




Within 
vv 1 LI III 1 


3 54E-06 


0 144 


c 


Sindex 


62 


ndpi 1 idpo 1 zot D 1 u j:?zu'+ 




H-JU:70H-DZ 




1 35929 


1 09E-05 


0 443 


c 


Sindex 


87 


l-l3nmai-iilQ7Afi-RTA-7fi1 Oft 
ndpi 1 la^JH^ / D 1 r\ / O 1 UD 








Within/R'^fiAQ 

vv ILI III 1/ OZOH-7 


1 97E-05 


0 801 


c 


4index 


52 


ARS-RFHI -NG*;-1 nfifll s 
/Anj Drvji— iiiuj 1 \juw 1 J 




61 1 99572 


RRM47 


Within 

V V 1 LI III 1 


2 03 E-06 


0 082 


c 


4index 


54 


ARS-RFfil -MGS-Rn54R 


7 


6434821 


r7H1Qnrf44 


Within 

V V 1 LI III 1 


2 79E-06 


0 113 


c 


4index 


42 


ttnj DrUL INUj 1 1 U30 1 


7 


1 5782979 


COL5A3 


Within 
vv ILI III 1 


5.50E-07 


2.23E-02 


c 


4index 


78 


Mnj DrVJL INUj 1 UHUjU 


7 


ft9R^Qc;Rn 

DZO Ji* JOU 


l_ J/ vf\ I r\ 1 


1 3832 


1 70E-05 


0 691 


c 


Sindex 


77 


BTB-0031 6650 


7 


6509S028 


GLRA1 


Within 


1 .67E-05 


0.679 


c 


Sindex 


63 


BTB-01 687547 


8 


20989026 


LOCI 01 905651 


328218 


1.13E-05 


0.458 


c 


Sindex 


48 


ARS-BFGL-NGS-9052 


8 


78009328 


FRMD3 


181107 


6.84E-06 


0.278 


c 


Sindex 


90 


ARS-BFGL-NGS-1 06379 


S 


113159018 


TSSCl 


90395 


2.24E-05 


0.909 


c 


Sindex 


52 


ARS-BFGL-NGS-1 5823 


9 


28887462 


PKIB 


Within 


8.1 5 E-06 


0.331 


c 


Sindex 


10 


BTB-01 33299S 


10 


73466092 


SLC38A6 


15901 


6.76E-07 


2.74E-02 


c 


Sindex 


29 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC38A6 


37359 


3.30E-06 


0.134 


c 


Sindex 


78 


ARS-BFGL-NGS-22113 


10 


73551579 


TMEM30B 


94651 


1 .68E-05 


0.6S3 


c 


Sindex 


26 


BTB-00471219 


11 


32709462 


NRXNl 


Within 


3.14E-06 


0.128 
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Trait 


Rank" 


SNP 


Chr. 


Positlon(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value* 


P value 
Bonferroni 


C 


Sindex 


44 


Hapmap25798-BTA-1 2638S 


11 


32731961 


NRXNl 


Witiiin 


6.07 E-06 


0.247 


C 


Sindex 


6 


Hapmap40257-BTA-91 91 6 


11 


32789048 


NRXNl 


Withiin 


3.25E-07 


1 .32E-02 


c 


Sindex 


93 


ARS-BFGL-NGS-S5007 


11 


53430229 


CTNNA2 


1292050 


2.37E-05 


0.964 


c 


Sindex 


72 


BTB-0107927S 


11 


5707S447 


REG3A 


436061 


1 .39E-05 


0.563 


c 


Sindex 


60 


BTB-01 079350 


11 


57107070 


REG3A 


4646S4 


9.51 E-06 


0.3S6 


c 


Sindex 


1 


Hapmap5 1 53 1 -BTA-98947 


11 


57522675 


REG3A 


8802S9 


2.08E-08 


8.44E-04 


c 


Sindex 


2 


BTB-01 328920 


11 


57639084 


LRRTM4 


792049 


4.34E-08 


1.76E-03 


c 


Sindex 


54 


ARS-BFGL-NGS-1 14087 


11 


64057850 


SPRED2 


293893 


8.21 E-06 


0.333 


c 


Sindex 


25 


ARS-BFGL-NGS-43985 


11 


64184454 


SPRED2 


420497 


3.00E-06 


0.122 


c 


Sindex 


40 


ARS-BFGL-NGS-91014 


11 


65493222 


a lH2orf66 


Witiiin 


4.88E-06 


0.198 


c 


6index 


8 


ARS-BFGL-NGS-1 1086S 


11 


75329413 


KLHL29 


Witiiin 


1.65E-05 


0.670 


c 


Sindex 


89 


UA-IFASA-2295 


11 


97754180 


RALGPSl 


Witiiin 


2.16E-05 


0.878 


c 


Sindex 


27 


BTA-119672-no-rs 


11 


102911946 


AK8 


Witiiin 


3.21 E-06 


0.130 


c 


Sindex 


71 


BTB-00490466 


12 


47836570 


DISS 


Witiiin 


1.37E-05 


0.557 


c 


4index 


81 


Hapmap42477-BTA-22799 


12 


48645720 


KLF12 


181742 


2.25E-05 


0.914 


c 


Sindex 


83 


Hapmap57649-rs2902241 4 


12 


6606S259 


GPC5 


2238S7 


1.85E-05 


0.753 


c 


6index 


2 


Hapmap381 87-BTA-l 05082 


12 


69216S40 


OCT 


286101 


7.38E-06 


0.300 


c 


Sindex 


4 


ARS-BFGL-BAC-137S8 


12 


69512332 


OCT 


Witiiin 


1 .79E-07 


7.27E-03 


c 


Sindex 


46 


ARS-BFGL-NGS-45730 


12 


70260457 


A8CC4 


49496 


6.28E-06 


0.255 


c 


4index 


55 


ARS-BFGL-NGS-1 3252 


12 


81184160 


GGACT 


25567 


3.04E-06 


0.123 


c 


6index 


6 


Hapmap53988-rs29024591 


13 


65855988 


EPB41L1/ACSS2 


23900/1 01 446S 


1.19E-05 


0.4S2 


c 


6index 


9 


Hapmap4071 2-BTA-33406 


13 


67101174 


BLCAP/ACSS2 


1353S/2259654 


1 .74E-05 


0.705 


c 


Sindex 


61 


ARS-BFGL-NGS-40278 


13 


67122063 


BLCAP/ACSS2 


Witllin/22S0543 


9.81 E-06 


0.398 


c 


6index 


3 


BPI-2 


13 


67833153 


L0CS14211/ACSS2 


Witllin/2991633 


7.51 E-06 


0.305 


c 


6index 


10 


ARS-BFGL-NGS-1 071 13 


13 


6795S189 


LOGS 1 4978/ ACSS2 


Witllin/31 16669 


1.83E-05 


0.743 


c 


eindex 


1 


Hapmap55254-rs2901 4939 


13 


69042143 


DHX35 


573959 


5.20E-06 


0.211 






41 


Hapmap30381 -BTC-005750 


14 


1463676 


C14H8orf33/DGAT1 


23690/331749 


5.62E-06 


0.228 




4index 


73 


ARS-BFGL-NGS-57820 


14 


1651311 


LOCW0294916/DGAT1 


Witllin/144114 


1.03E-05 


0.419 




4index 


75 


ARS-BFGL-NGS-4939 


14 


1801 116 


DGAT1 


Witiiin 


1.15E-05 


0.469 




Sindex 


88 


ARS-BFGL-NGS-4939 


14 


1801 116 


DCAT1 


Witiiin 


2.00E-05 


0.812 




4index 


47 


Hapmap52798-ss46526455 


14 


1923292 


MAF1/DGAT7 


Witliin/1 1S454 


1 .08E-06 


4.39E-02 


Q 


4mdex 


80 


Hapmap30986-BTC -056068 


14 


10346734 


EFR3A 


64293 


1.97E-05 


0.799 




Sinclcx 


51 


BTB-00557532 


14 


24643266 


XKR4 


28554 


7.35E-06 


0.299 




Sindsx 


82 


ARS-BFGL-BAC-1 0245 


14 


31819743 


PDE7A 


Witiiin 


1 .80E-05 


0.731 




Sindex 


34 


ARS-BFGL-NGS-2025 


14 


65452336 


ZNF706 


103905 


3.79E-06 


0.154 




Sindex 


8 


Hapmap35102- 
BES3_Contig324_378 


14 


66523007 


RGS22 


Witiiin 


5.63E-07 


2.29E-02 


c 


Sindex 


58 


UA-IFASA-4785 


14 


71096693 


C14H8orf37 


300038 


9.36E-06 


0.3S0 


c 


Sindex 


79 


BTB-0129621S 


15 


17344116 


ALKBH8 


21050 


1.77E-05 


0.720 


c 


Sindex 


39 


ARS-RFGI -MGS-^^7n4 


1 5 


1 

1 Oj J J 


l_ / Dn 1 1 Ul lOD 


Within 

vv III III 1 




0 1 96 


c 


Sindex 


37 


ARS-BFGL-NGS-1 2066 


15 


27291268 


BUD13 


542743 


4.33E-06 


0.176 


c 


Sindex 


43 


BTA-36518-no-rs 


15 


32983903 


SORLl 


228564 


6.04E-06 


0.245 


c 


Sindex 


81 


Hapmap56991 -rs2901 00S3 


15 


81595546 


LOCS38839 


52784 


1 .78E-05 


0.724 


c 


Sindex 


68 


Hapmap46697-BTA-381 71 


16 


2899256 


NUAK2 


Witiiin 


1 .24E-05 


0.501 


c 


Sindex 


7 


ARS-BFGL-NGS-6S533 


16 


60615012 


LOG 101 902340 


3599S 


S.39E-07 


2.19E-02 


c 


6index 


7 


BTB-01 090859 


16 


6324S646 


XPRl 


51157 


1.51 E-05 


0.615 


c 


Sindex 


64 


ARS-BFGL-NGS-36880 


16 


73736551 


SLC30A1 


28425 


1.14E-05 


0.465 


c 


4index 


83 


ARS-BFGL-NGS-38696 


17 


1839S611 


MGST2 


29230 


2.29E-05 


0.929 
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Trait 


Rank^ 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value^ 


P_ value 
Bonferroni 


— 

CI 4inclGX 


64 


/An J DrUL l\U J 00?Z J 


1 7 


1 RQ1 flQIA 
1 1 1 H 


CCRN4L 


Within 

vv III III 1 


5 55E-06 


0 225 


CI 4index 


50 


ADC RP("I Mf~C imQ^? 
rtnb-bruL-lNu J- 1 UzyjJ 


1 7 


1 1 CTQ/I 

1 yj 1 jzy4 


LLn/V4L 


ACiA ^AQ 
4U4.50O 


1 .81 E-06 


0.074 


CI 4index 


39 


D 1 A-y 1 D / D-no-rs 


1 7 


zuuyyoj/ 




291 555 


2.96E-07 


1 .20E-02 


CI 4index 


79 


napmapD 1^4 .d-d i a-'Hjd i y 


1 7 


zuyyoo4/ 




380577 


1 .72E-05 


0.698 


CI 4indGX 


53 


D 1 D U 1 JOJZ.Ui? 


1 7 


1 U 1 / O 




53460 


2 79 E-06 


0 113 


1 nH fiv 


35 




1 7 


449631 24 




514 


U ■ " \J ^ VJ o 


2.62E-03 


k_ 1 Oil lUcA 


1 9 


/Anj DrVJL INVJJ / 1 1 ID 


1 7 


ftRnn'^'^.in 

DOUUZ JH-U 


IVl 1 \J 1 OD 


Within 

vv III III 1 


2 41 E-06 


0 098 


L_ 1 0 1 1 1 U C A 


36 


ARc;-RF(^l 


1 7 


ftRilQflAc;^ 


TPST2 


Within 
vv 1 11 MM 


A ^flF-flfi 


0 175 


f"1 RinHov 
^ 1 Oil ILltrA 


74 


RTR-ni TOnHdfi 

D 1 D u 1 / ?UOM-0 


20 


1 32631 57 


SREK1 


1 28568 


1 44E-05 


0 586 


CI Rinrlp>v 
^ 1 Oil lUcA 


85 


AR^i-RFf^l -RAr-'JdfiQ 


20 


J jM-j J 1 QU 


HFA TR 7ft :? /nycTi 


'>l^n7 /^RAA^A 

i.£.J\J 1 1 JOH'+jH- 


1 93E-05 


0 784 


^ 1 Oil lUUA 


56 


AR<;-RF(^I -Nlf^<^-7fi7^fi 

/An J DrVJL 1 \U J / D/ JO 


20 


33491 273 


HF4 Tf? 7ft ? /nycTi 


Within/fi/l'?'^jl7 

vv III III 1/ OM-Z DM-/ 


9 20E-06 


0 373 


CI RinHov 
L_ 1 0 1 1 1 U t: A 


49 


D 1 D U 1 HZ jD J J 


20 


DO J / OZUU 


jrcrz./ rni.n 


WithinMQ'^flilft 

vv 1 11 III l/H!? JUM-D 


D.^DC UO 


0 282 


L 1 omuex 


42 


D 1 D-U I4ZDD/D 


20 


DtSoUODDD 


CDPP'J/DDI D 
jrtrz/rnLn 


vviinin/4DDoy.5 


6.02 E-06 


0.244 


n RinHpY 

^_ 1 0 1 1 lU c:A 


20 


H3nm3n3n57n-RTA-1 S?77R 

nuLliiialJ3VJJ/L/ Dl/A IJZ//0 


20 


38761 71 1 

30 / L; 1 / 1 1 


3t L.i Z.I rtlLft 


1 54961 /31 1 

1 _J*T 1 / J 1 1 J J J 


2.55E-06 


0 1 04 


n RinHpY 

V_ 1 01 1 lU ^A 


80 


AR*;-RF(^I -Nf^*;-QQ71 fi 


21 


63560939 


VRK1 


216337 


1 78E-05 


0 723 


n RinHov 

1 Oil lUcA 


53 


ARC;-RF(^I -Mr;'s-^Q4'^Q 
/An J DrVJL InU J j!7H-Jj 


22 




/ v/nZ LyZ 


Within 

vv 1 11 III 1 


8 1 5 E-06 


0 331 


L 1 oinaex 


69 


ndpmdpD^DDo-rszyuuyDyis 


24 


zy 1 0/0U4 




vviinin 


1 .27E-05 


0.51 7 


CI 6indsx 


4 


AR*^-RF(^I -Mr;'s-4'^ft7Q 

/Anj DrVJL InU J M-JD/ 


24 


H-Z JOZ DU J 


APCDD 1 


Within 
vv III III 1 


7 82 E-06 


0 318 


CI Gindcx 




rid pi 1 idpo 1 ZDU D 1 u i j jz / 


25 


2224930 


ZG16B 


Within 
vv III III 1 


8 51 E-06 


0 346 


CI 4indcx 


51 


rid pi 1 (dp jh-uom- 1 bz?u i i i?30 


26 


5526925 


PCDH15 


Within 
vv III III 1 


2 02 E-06 


0 082 


CI 4indGX 


77 


RTR-fll n77Q^Q 

D 1 D U 1 U/ / 


26 


7ftRc;i 1 n 

/ DO J 1 1 U 


PRKG1 


Within 
vv III III 1 


1 35E-05 


0 548 


CI 4indGX 


27 


AR<;-RFril -Nr;'s-1^74ft 

/An J DrUL 1 \VJ J 1 J / H-O 


26 


QRfiftQAfl 

^ODO!7H-U 


RNLS 


651 37 


3.71 E-09 


1 .51 E-04 


C 1 AinH ov 
^ 1 T-l 1 lU trA 


72 


ndlJI 1 ICIIJJO 1 OJ 1 DZ.^KJZ-Z.Z-^'-r 


26 


10009077 


RNLS 


Within 

V V 1 LI III 1 


1 02E-05 


0 415 


C 1 Ai nH ov 
v_ 1 *tI 1 iU trA 


35 


RTA-filQ^I-nn-rt; 
Di/Aui_7zi iiur^ 


26 


1 L/Z-J JZ JO 


LIPJ 


7406 


6.99E-08 


2.84E-03 


CI 4indGX 


20 


AR<;-RF(^I -HC-i^-lMQA 

/An J DrUL INU J Z 1 / 


26 


1 KJD^ 1 30Z 


\-\JK- 1 \J\Jdj\jDD / / i-Irn 


1 D jUi?/ DDZOO 


4.20E-1 0 


1 .70E-05 


C 1 4indsx 


68 


rid pi 1 idp J J J 1 bz^u 1 Doou 


26 


1 UDOi? J / ^ 


ACTA2 


9731 


O. jOC UO 


0 339 


CI 4indGX 


76 


RTA-1 1 1 R^7-nn-t-c 
D 1 r\ 1 1 IOj/ IIL; lb 


26 


1 UO 1 JDOO 


FAS 


67560 


1 17E-05 


0 476 


CI 4indsx 


56 


D 1 D UU^ZM-U 1 J 


26 


\ U^ZZUO 1 


CH25H 


541 32 


3 05 E-06 


0 1 24 


CI 4indGX 


48 


RTA-nR77c;-rc9Qn99^^'5 
D 1 r\ UO/ / J 1 r>Z7UZZ3 jZ 


26 


1 1 9ni 1 QR 

1 1 ZU 1 1 70 


-JH_ / Or\ iZ 


39720 


1 .09E-06 


4.42E-02 


CI 4indGX 


57 


ARci-RFni -Mr;<;-'^^i i 

/Anj DrVJL INUj J J 1 1 J 


26 


1 1 528933 


KIF20B 


61 222 


3 1 0E-06 


0 1 26 


CI 4indGX 


25 


AR<;-RF(^I -Mr;<^-fi^R^^ 

/An J DrUL INU J OjO J J 


26 


1 1 3HZODO 




67662 


2.50E-09 


1 .01 E-04 


C 1 Ai nH ov 
v_ J *tI 1 lu trA 


61 


AR*^-RFfil -Mr;'s-19381 

/Anj DrVJL INUj 1 ZJO 1 


26 


1 :?:?00Q48 

1 ZZVjU_7*tO 


1 onnnidi 74? 


24477 


4 76E-06 


0 1 93 


C 1 4.inH ov 
V. 1 T^l 1 lU trA 


1 7 


RTR-niQnR417 


26 


1 ?:?6R4!?7 

1 ZZL;OH-Z / 




15486 


1 .71 E-1 0 


6.94E-06 


C 1 AinH ov 
V_ 1 M-| 1 lu trA 


67 


RTR-fll 841 68'? 

D 1 D U 1 OH* 1 UOZ 


26 


12295284 




42343 


8 27E-06 


0 336 


CI 4index 


84 


ADC RP("i Mf~c A^^A^^ 


26 


1 z3o4 1 1 y 


nln/ 


1 0881 


2.32 E-0 5 


0.942 


CI 4indGX 


23 


rid pi 1 idpDZO 1 / 1 bz?u i i i'oy 


26 


141 55229 


HHEX 


291 60 


1 .34E-09 


5.45E-05 


CI 4indGX 


69 


AR*;-RF(^I -Mf^'^-R'^RfiA 

/An J DrvjL INU J O j0D'+ 


26 


14532797 


CYP26A 1 


69001 


R ^RF-flfi 

O. jOC UO 


0 340 


CI 4!ndGX 


74 


ARCl-RF^il -Mri's-1 1047^^ 

/Anj DrvjL InU J 1 1 UH-/ J 


26 


1 t^ftflAfi^l 

1 JDUM-O J 1 


PLCEl 


Within 
vv III MM 


1 04 E-0 5 


0 421 


Cl AinHfiv 

V. 1 HII lUCA 


58 


RTR-fin7nfiR^R 

D 1 D UU/ UDOjO 


26 


1 JOZH- 1 M- 1 


TBC1D12 


Within 

vv III III 1 


3 31 E-06 


0 1 34 


C14indGX 


70 


BTB-00927439 


26 


16315378 


CYP2C19 


20815 


9.08E-06 


0.369 


CMindGX 


18 


ARS-BFGL-NGS-29299 


26 


16614068 


PDUM1/S0RBSJ 


Within/73061 


2.06E-10 


8.36E-06 


C14indGX 


49 


Haprrap41595-BTA-60800 


26 


16791783 


SORBS 7 


Within 


1 .09E-06 


4.43E-02 


CMindGX 


28 


Hapmap58930-rs2901 0490 


26 


16822073 


SORBS 7 


Within 


5.90E-09 


2.40E-04 


C14indGX 


59 


ARS-BFGL-NGS-1 06959 


26 


17225652 


CC2D2B 


Within 


3.80E-06 


0.154 


CMindGX 


63 


ARS-BFGL-NGS-1 13660 


26 


1 7246984 


CC2D2B 


Within 


5.20E-06 


0.211 


CMindGX 


71 


ARS-BFGL-NGS-25217 


26 


17307507 


CCNJ 


13017 


9.96E-06 


0.404 


CMindex 


66 


ARS-BFGL-NGS-1 14539 


26 


18808408 


SFRP5 


21032 


7.36E-06 


0.299 


CMindGX 


82 


ARS-BFGL-NGS-97471 


26 


18882047 


CRTACI 


109 


2.27E-05 


0.921 
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Table 5. Cont. 





Trait 


Rank" 


SNP 


Chr. 


Positlon(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value" 


P_value 
Bonferroni 


C14index 


34 


ARS-BFGL-NGS-41056 


26 


18906121 


CRTACl 


Withiin 


4.78E-08 


1 .94E-03 


C14index 


41 


ARS-BFGL-NGS-1 16902 


26 


18967997 


CRTACl 


Within 


4.05E-07 


1 .64E-02 


C14index 


43 


ARS-BFGL-NGS-25126 


26 


18994785 


CRTACl 


Within 


S.96E-07 


2.42E-02 


C14index 


62 


ARS-BFGL-NGS-71584 


26 


20290497 


GOTl/SCDl 


Within/842247 


5.1 1 E-06 


0.207 


CMIndex 


M 


ARS-BFGL-NGS-23064 


26 


20365711 


NKX2-3/SCD1 


34516/767033 


2.95E-1 1 


1 .20E-06 


CMindex 


15 


ARS-BFGL-NGS-77668 


26 


20393457 


NKX2-3/SCD1 


6770/74448S 


6.02E-1 1 


2.44E-06 


CMindex 


33 


ARS-BFGL-NGS-2464 


26 


20444634 


SLC25A28/SCD1 


21194/688110 


3.9eE-08 


1.61 E-03 


CMindex 


19 


BTB-0a930925 


26 


20474308 


SLC25A28/SCD1 


Within/65S436 


2.61 E-10 


1.06E-05 


CMindex 


24 


ARS-BFGL-NGS-39397 


26 


20716721 


DNMBP/SCDJ 


Within/416023 


1 .54E-09 


6.25E-05 


CMindex 


22 


BTB-00930720 


26 


20903573 


LOC5U498/SCD1 


Within/244744 


6.89E-10 


2.80E-05 


CMindex 


31 


Hapmap4641 l-BTA-1 5820 


26 


20984335 


CHUK/SCDl 


Within/1 4S409 


1.10E-08 


4.45E-04 


CMindex 


16 


Hapmap31825-BTA-1 58647 


26 


21056547 


PKD2U/SCD1 


Within/76197 


1.56E-10 


6.35E-0e 


CMindex 


10 


Hapmap33073-BTA-1 62864 


26 


21180893 


SCD1 


32576 


3.04E-1 2 


1.23E-07 


CMindex 


1 


BTB-00931481 


26 


21226405 


WNT8B/SCD1 


14100/7808S 


6.91 E-1 7 


2.80E-12 


CMindex 


13 


ARS-BFGL-NGS-1 10077 


26 


21322557 


HinAN/SCDl 


22399/174240 


2.04E-1 1 


8.28E-07 


CMindex 


11 


ARS-BFGL-NGS-1 08305 


26 


21363670 


SCD1/HIFIAN 


63512/215353 


6.46E-12 


2.62E-07 


CMindex 


5 


BTB-00931586 


26 


21409429 


PAX2/SCD1 


61334/261112 


3.39E-13 


1 .38E-08 


CMindex 


4 


ARS-BFGL-NGS-1 14149 


26 


21702714 


LZTS2/SCD7 


656/564769 


3.10E-13 


1 .26E-08 


CMindex 


9 


ARS-BFGL-NGS-1 16481 


26 


21977581 


LOC100847491/SCD1 


10062/829264 


1.20E-12 


4.87E-08 


CMindex 


7 


Hapmap24832-BTA-1 38805 


26 


22016380 


BTRC/SCDl 


Within/86S063 


6.31 E-1 3 


2.56E-08 


CMindex 


6 


ARS-BFGL-NGS-6259 


26 


22059103 


BTRC/SCD1 


Within/910786 


S.97E-13 


2.42E-08 


CMindex 


8 


BTB-00932332 


26 


2211S554 


BTRC/SCDl 


Within/970237 


7.55E-13 


3.06E-08 


CMindex 


2 


ARS-BFGL-NGS-1 07403 


26 


22889812 


NFKB2 


1586 


2.62E-15 


1.06E-10 


CMindex 


21 


Hapmap48222-BTA-1 22240 


26 


23641881 


C26HWorf26 


Within 


5.68E-10 


2.31 E-05 


CMindex 


30 


Hapmap49372-BTA-91 009 


26 


23689229 


C26HWorf26 


3017 


6.47E-09 


2.63E-04 


CMindex 


32 


BTA-60918-no-rs 


26 


23853334 


CNNM2 


Within 


3.34E-08 


1.36E-03 


CMindex 


12 


BTA-60935-no-rs 


26 


23876476 


CNNM2 


Within 


1 .98E-1 1 


8.03E-07 


CMindex 


40 


ARS-BFGL-NGS-2180 


26 


24477962 


SH3PXD2A 


Within 


3.20E-07 


1.30E-02 


CMindex 


29 


ARS-BFGL-NGS-1 092 


26 


24531763 


SH3PXD2A 


Within 


6.21 E-09 


2.52E-04 


CMindex 


65 


ARS-BFGL-NGS-1 81 94 


26 


24575207 


SH3PXD2A 


Within 


5.58E-06 


0.227 


CMindex 


3 


ARS-BFGL-NGS-1 181 89 


26 


24786731 


SLK 


Within 


6.55E-14 


2.66E-09 


CMindex 


46 


UA-IFASA-4715 


26 


25314352 


CCDC147 


27285 


9.22E-07 


3.74E-02 


CMindex 


44 


BTB-00935537 


26 


26585557 


S0RCS3/ECHS1 


187266/719949 


7.48E-07 


3.04E-02 


CMindex 


37 


Hapmap28763-BTA-1 6232S 


26 


26757136 


S0RCS3/ECHS1 


358845/891528 


7.79E-08 


3.16E-03 


CMindex 


60 


BTA-S7355-no-rs 


26 


27251857 


SORCSl 


558341 


4.71 E-06 


0.191 


CMindex 


45 


ARS-BFGL-NGS-1 448 


27 


37357125 


HOOK3/AGPAT6 


Within/1128138 


7.66E-07 


3.11E-02 


CISindex 


11 


ARS-BFGL-NGS-1 10992 


2S 


20421361 


REEP3 


614577 


9.76E-07 


3.96E-02 


ClSindex 


18 


ARS-BFGL-NGS-4865 


2S 


28011033 


CDH23 


Within 


2.38E-06 


0.097 


CISindex 


86 


ARS-BFGL-NGS-1 2970 


29 


40646639 


PPPm32/FADS1 


2220/292226 


1.94E-05 


0.7S8 


CISindex 


57 


ARS-BFGL-NGS-1 1681 


X 


3622615 


SLC25A43 


24203 


9.32E-06 


0.378 


CISindex 


67 


Hapmap48540-BTA-97806 


X 


8700029 


ODZl 


Within 


1.23E-05 


0.498 


CISindex 


70 


Hapmap50046-BTA-58882 


X 


82022276 


MIR374B 


1635 


1.33E-05 


0.540 


CISindex 


17 


Hapmap60551-rs29017241 


X 


82281306 


XIST 


Within 


2.17E-06 


0.0S8 


CISindex 


38 


Hapmap60664-rs2901 7374 


X 


107043386 


CASK 


1 1 38S9 


4.61 E-06 


0.1 S7 


CISindex 


66 


Hapmap49563-BTA-30596 


X 


120716539 


lURAPU 


Within 


1.21 E-05 


0.490 



Note: see note to Table 2. 

doi:1 0.1 371 /journal.pone.00961 86.t005 
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Table 6. Genome-wise and suggestive significant SNPs for sum of fatty acid traits. 





Trait 


Rank" 


SNP 


Chr. 


Position(bp) 


Candidate Gene 


Distance(bp) 


P_value' 


Bonferroni 


SFA 


13 


Hapmap42233-BTA-49670 


1 


82559884 


CmSorfJO/EHHADH 


14106/38733 


6.58E-06 


0.267 


SPA/UFA 


14 


Hapmap42233-BTA-49670 


1 


82559884 


CmSorfJO/EHHADH 


14106/38733 


4.96E-06 


0.202 


UFA 


33 


ARS-BFGL-NGS-13938 


2 


111650859 


M0GAT1/ACSL3 


19295/146311 


2.25E-05 


0.912 


SFA 


8 


Hapmap42304-BTA-73062 


5 


5049476 


KRR1/OSBPL8 


14567/770655 


5.51 E-06 


0.224 


SFA/UFA 


12 


Hapmap42304-BTA-73062 


5 


5049476 


KRR1/OSBPL8 


14567/770655 


4.17E-06 


0.169 


UFA 


28 


Hapmap42304-BTA-73062 


5 


5049476 


KRR1/OSBPL8 


14567/770655 


2.02E-05 


0.819 


SFA/UFA 


32 


ARS-BFGL-NGS-38038 


5 


27992179 


NR4A1 


Within 


1 .80E-05 


0.731 


UFA 


26 


ARS-USMARC-624 


5 


28859701 


CSRNP2 


Within 


1 .87E-05 


0.760 


SFA 


10 


ARS-BFGL-NGS-8796 


5 


29095603 


LOC510716 


Within 


6.04E-06 


0.245 


SFA/UFA 


8 


ARS-BFGL-NGS-8796 


5 


29095603 


LOC510716 


Within 


3.04E-06 


0.123 


UFA 


3 


ARS-BFGL-NGS-8796 


5 


29095603 


LOC510716 


Within 


1 .44E-06 


0.058 


SFA/UFA 


22 


ARS-BFGL-NGS-35179 


5 


34325053 


SCAF1 1 


62336 


1.21 E-05 


0.493 


UFA 


20 


ARS-BFGL-NGS-35179 


5 


34325053 


SCAFl 1 


62336 


1.43E-05 


0.579 


SFA 


17 


ARS-BFGL-NGS-69056 


5 


42285835 


CPNE8 


135184 


1.11 E-05 


0.449 


SFA/UFA 


29 


ARS-BFGL-NGS-69056 


5 


42285835 


CPNE8 


135184 


1.62 E-05 


0.657 


UFA 


8 


ARS-BFGL-NGS-69056 


5 


42285835 


CPNE8 


135184 


4.25E-06 


0.173 


UFA 


27 


Hapmap52463-r52902583 1 


5 


45473334 


LOCI 01 905276 


65002 


1.93 E-05 


0.784 


SFA 


11 


Hapmap39862-BTA-74478 


5 


85672503 


BCATl 


74712 


6.22E-06 


0.253 


SFA/UFA 


6 


Hapmap39862-BTA-74478 


5 


85672503 


BCATl 


74712 


2.26E-06 


0.092 


UFA 


4 


Hapmap39862-BTA-74478 


5 


85672503 


SCAT! 


74712 


1.92 E-06 


0.078 


SFA/UFA 


16 


ARS-BFGL-NGS-99256 


5 


104714350 


VWF/OLR7 


Within/4458707 


5.47E-06 


0.222 


UFA 


12 


ARS-BFGL-NGS-99256 


5 


104714350 


VWF/OLR7 


Within/4458707 


6.47E-06 


0.263 


UFA 


16 


BTB-00246150 


6 


20993424 


PPA2 


Within 


9.44E-06 


0.383 


SFA/UFA 


34 


Hapmap26001 -BTC-03881 3 


6 


44926243 


PPARCC1A 


Within 


1 .97E-05 


0.799 


UFA 


25 


Hapmap26001 -BTC-03881 3 


6 


44926243 


PPARCC1A 


Within 


1.85 E-05 


0.750 


UFA 


31 


BTB-00316291 


7 


64892251 


SPARC 


Within 


2.18E-05 


0.887 


SFA 


15 


BTB-003 16348 


7 


64939808 


ATOXl 


5940 


8.62E-06 


0.350 


SFA/UFA 


17 


BTB-003 16348 


7 


64939808 


ATOXl 


5940 


6.58E-06 


0.267 


UFA 


15 


BTB-003 16348 


7 


64939808 


ATOXl 


5940 


9.25E-06 


0.375 


SFA 


9 


BTB-003 16650 


7 


65098028 


GLRAl 


Within 


6.02E-06 


0.245 


SFA/UFA 


18 


BTB-003 16650 


7 


65098028 


GLRAl 


Within 


7.85E-06 


0.319 


UFA 


17 


BTB-003 16650 


7 


65098028 


GLRAl 


Within 


9.98E-06 


0.405 


UFA 


30 


BTB-02040446 


8 


30320626 


NFIB 


164574 


2.16E-05 


0.878 


UFA 


22 


HapmapSOl 26-BTA-83733 


9 


55449737 


LOCW1907134 


18763 


1.73E-05 


0.702 


SFA 


24 


BTB-0 18665 13 


10 


70231445 


SLC35F4 


Within 


1.73E-05 


0.701 


SFA 


29 


BTB-01203179 


10 


72694329 


DHRS7 


53135 


2.03E-05 


0.825 


SFA/UFA 


28 


BTB-01203179 


10 


72694329 


DHRS7 


53135 


1.61 E-05 


0.653 


UFA 


10 


BTB-01203179 


10 


72694329 


DHRS7 


53135 


5.18E-06 


0.210 


SFA 


2 


BTB-01 332998 


10 


73466092 


SLC38A6 


15901 


8.08E-08 


3.28E-03 


SFA/UFA 


2 


BTB-01 332998 


10 


73466092 


SLC38A6 


15901 


3.54E-07 


1 .44E-02 


UFA 


2 


BTB-01 332998 


10 


73466092 


SLC38A6 


15901 


4.87E-07 


1 .98E-02 


SFA 


1 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC38A6 


37359 


6.07E-08 


2.46E-03 


SFA/UFA 


1 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC38A6 ___ 


37359 


2.64E-07 


1 .07E-02 


UFA 


1 


ARS-BFGL-NGS-4783 


10 


73487550 


SLC38A6 


37359 


4.01 E-07 


1.e3E-02 


SFA 


7 


BTB-01 501 723 


10 


73655817 


TMEM30B 


6161 


5.31 E-06 


0.216 


SFA/UFA 


37 


BTB-01 501 723 


10 


73655817 


TMEM30B 


6161 


2.28E-05 


0.924 


UFA 


34 


BTB-01 501 723 


10 


73655817 


TMEM30B 


6161 


2.25E-05 


0.912 


SFA 


19 


BTB-01 079278 


11 


57078447 


REG3A 


435723 


1 .39E-05 


0.565 


SFA/UFA 


36 


BTB-01 079278 


11 


57078447 


REG3A 


435723 


2.12E-05 


0.859 
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Table 6. Cont. 















Nearest Gene/ 




Raw 


P_ value 


Trait 


Rank'' 


SNP 


Chr. 


Position(bp) 


Candidate Gene 


Distance(bp) 


P_value^ 


Bonferroni 


SFA 


23 




1 1 


57107070 




454346 


1 66E-05 


0 675 


SFA 


3 


RTA-1 1 Qfi77-nn-rt; 

DIM 1 J / Z. TIL/ 1 ^ 


1 1 


1 0991 1 946 

1 UZ_? 1 1 _?*tLJ 


AK8 


Within 

VV 1 LI III 1 


1 34E-06 


0 054 


cpA/l )FA 


10 


RTA-1 1 Qfi7!?-nn-r<; 

D 1 rA II / Z. 1 1 U 13 


1 1 


1 0791 1 946 

1 UZ_7 1 1 _7*tU 


AK8 


Within 

V V 1 LI III 1 


3 77E-06 


0 1 53 


UFA 




RTA-1 1 QfiT'J-nn-rc 

D 1 tt 1 1 ^\J/ £. \ ID 1 b 


1 1 


1 n9Q1 1 QAf, 
\ uzy 1 1 ^^+0 


AK8 


Within 

VVI Llllll 


M-.UjC UD 


0 164 


UFA 


1 8 


D 1 D U 1 Z jO^Ui* 


1 2 


502951 1 2 


TBC1D4 


382864 


1 1 7E-05 


0 476 


UFA 


38 


RTR-ni QRn4H7 

D 1 D U 1 i?OUHOZ 


1 2 


JUM-J I Z03 


TBC1D4 


226687 


2 38E-05 


0 964 




31 


MnJ DrVJL InVJ J / UZUD 


1 3 


M-OOZZD J J 


FFRMT1 


Within 

VVI LI III 1 


1 65E-05 


0 670 


UFA 


37 


AR^;-RFr;i -Mr;^-7n7nfi 


1 3 


'+ODZZD J J 


FFRMT1 
rcr\ivi 1 1 


Within 

VVI LI III 1 


2 31 E-05 


0 937 


cpA/l IFA 
J r r\/ u rrt 


25 


AR*i-RF(^l -M(^^-'^7R7n 


1 4 


1 651 31 1 


L\J\^ 1 \j\Jz.yHy 1 o/ lyxj/i i i 


Within/1 AA^ 1 A 

VVI LI III 1/ 1 H-H- 1 1 


1 41 E-05 


0 572 


SFA 


5 


ARS-RFf^i -Mf^S-4Q^q 


14 


1 ROI 1 1 6 

1 0\J 1 1 1 D 


DGAT1 


Within 

VV 1 LI III 1 


T RnF-n6 

Z.OUC UD 


0 114 


SFA/l IFA 


3 


ARS-RFf;i -Mf;<i-4Q3q 


14 


1 ROI 1 1 6 

1 0\J 1 1 1 D 


DGA TI 


Within 

VV 1 LI III 1 


1 .15E-06 


4.65E-02 


UFA 




ARS-RFf^l -Mr;'^-4Q3q 


14 


1 Qfll 1 1 6 
1 0\J 1 1 1 D 


DGA TI 


Within 

VV 1 LI III 1 


4 78E-06 


0 1 94 


cpA/l IFA 
J r r\/ u r r\ 


21 


AR'^-RFril -Mf^Ci-l n7^7Q 

rtn J DrUL InO J 1 U / J / y 


1 4 


ZU JH-H- J / 


L\-/i^/ OOy OO/ L/tJri I I 


'+DU/ ZH^D 1 y 


1 1 7E-05 


0 475 


cpA/I IFA 


1 3 


AR^i-RFr;! -Mri'^-i fifiiiRn 


14 


'+JDH-7 JZ 


TRAPPC9 


Within 

VVI LI III 1 


H-.DOC UD 


0 1 90 


UFA 


14 


AR^^-RFrii -Mr;<^-i nn^fin 


1 4 


A'^f^AQ^l 

'+JDH-3 JZ 


TRAPPC9 


Within 
VVI LI nil 


7 52E-06 


0 305 


SFA 


16 


AR^i-RFr;! -Mr;<^-1 1 ^9Q^ 


1 4 


7797jl^Rfi 


SLC2A5 


3495 


1 03 E-05 


0 419 


cpA/I IFA 
J r r\/ u r r\ 




AR'^-RFril -Mf^<^-1 1 ^^Q^ 


1 4 


7777il^Rft 

/ / Z/ HjOD 


SLC2A5 


3495 


1 41 E-06 


0 057 


UFA 


19 


AR'^-RFr;i -Mri'^-i i ^^q^ 


1 4 


7797il^Rfi 

/ / Z/ HjOD 


SLC2A5 


3495 


1 36E-05 


0 551 


SFA 


22 


ARS-RFr;i -Mr;<i-i s4ri 


1 8 


5661 1 355 




Within/RR9940 

V V 1 LI III 1/ OO _7 -^T-U 


1 55 E-05 


0 630 


SFA/l IFA 


20 


ARS-RFf^l 54R1 


1 8 


5661 1 355 

_J\JD I 1 


AP7Al/<iPHK7 


Within/RR9940 
VV 1 LI III 1/ \jQyy-v\j 


9 71 E-06 


0 394 


UFA 


24 


ARS-RFr;i -Mr;<i-i s4ri 


1 8 


5661 1 355 

_JLjD 1 1 


AP?A U*iPHK7 


Within/RR9940 
V V 1 LI III 1/ GQyy-v\j 


1 79E-05 


0 725 


SFA 


20 


AR*i-RFf;i -Nf;9-7^^Qn 


1 9 


M-ZOD/ OJ 


yujj 1 


6533 


1 44E-05 


0 586 


cpA/l IFA 
J r r\/ u rrt 


1 9 


AR^i-RFr;! -M(^'^-7R^Qn 


1 9 


M-ZOO/ OJ 




6533 


y.yjH C UD 


0 391 


SFA 




AR's-RFr;i -Mr;<^-R7^fiR 


1 9 


/ / DZOZU 


riQNI 7nrf/^7 
i_ ( yn 1 /ui lo/ 


67914 


3 15 E-06 


0 1 28 


cpA/l IFA 


g 


AR'^-RFr;i -Mr;<^-R7^fiR 


1 9 


/ / DZOZU 


i_ ( yn 1 / UI lo/ 


67914 


3 72 E-06 


0 151 


UFA 


1 3 


AR*^-RF(^I -M(^^-R7^fiR 


1 9 


/ / DZOZU 


V. 1 y" 1 / UI lo/ 


6791 4 


D.DH C UD 


0 270 


SFA 


21 


AR*^-RF^;l -Mr;<^-i i r^^q 


20 


33471 38 


v_Z Un DUI 1 DU 


J juouh- 


1 54E-05 


0 627 


SFA/UFA 


15 


ARS-BFGL-NGS-1 18339 


20 


3347138 


C20H5orf50 


1 50804 


5.31 E-06 


0.215 


UFA 


23 


ARS-BFGL-NGS-1 18339 


20 


3347138 


C20H5orf50 


1 50804 


1.75 E-05 


0.709 


SFA/UFA 


24 


BTA-50482-no-rs 


20 


36336225 


EGFLAM 


157753 


1.35E-05 


0.548 


UFA 


36 


BTA-50482-no-rs 


20 


36336225 


EGFLAM 


157753 


2.31 E-05 


0.936 


SFA 


14 


ARS-BFGL-NGS-1 16806 


20 


36450009 


GDNF 


1 80881 


7.68E-06 


0.312 


cpA/l IFA 


1 1 


AR'^-RFril -M(^<^-1 1 fiRflfi 


20 




GDNF 


1 anRH 1 

1 OUOO 1 


3 92 E-06 


0 159 


UFA 


1 1 


AR'^-RFr;i -Mri'^-i i fiROA 


20 


JOH-JUUU3 


GDNF 


1 OUOO 1 


5 23 E-06 


0 212 


SFA 


30 


Mn J U JlVlrtrW- rdl fcrl 1 L 

DQ990835-rs29012811 


20 


jO J / U JZy 


GDNF 


60361 


2 31 E-05 


0 939 


SFA/UFA 


27 


ARS-USMARC-Parent- 
DQ990835-rs29012811 


20 


36570529 


GDNF 


60361 


1 .56E-05 


0.632 


UFA 


35 


ARS-USMARC-Parent- 


20 


36570529 


GDNF 


60361 


2.29E-05 


0.931 


SFA 


27 


ARS-RFGI -MG^-17fi7fi 


20 


3901 79R5 


PRLR 


55261 


1 95 E-05 


0 790 


SFA/l IFA 


23 


ARS-RFGI -MG<i-17fi7fi 


20 


3901 79R5 


PRLR 


55261 


1 29E-05 


0 523 


UFA 


39 


ARS-BFGL-NGS-1 7676 


20 


3901 7985 


PRLR 


55261 


2.45 E-05 


0.994 


SFA/UFA 


38 


ARS-BFGL-NGS-55739 


20 


39787788 


OQTNFS 


Within 


2.38E-05 


0.965 


SFA 


25 


BTB-00783271 


20 


41201777 


SUB) 


17036 


1.73 E-05 


0.702 


SFA/UFA 


26 


BTB-00783271 


20 


41201777 


SUB1 


17036 


1 .54E-05 


0.626 


UFA 


29 


BTB-00783271 


20 


41201777 


SUB1 


17036 


2.09E-05 


0.850 


SFA 


12 


BTB-01 583562 


20 


55425112 


LOO 01 905359 


Within 


6.28E-06 


0.255 


SFA 


4 


Hapmap44836-BTA-51861 


21 


20990602 


ABHD2/PLIN1 


8117/512222 


1 .76E-06 


0.071 


SFA/UFA 


4 


Hapmap44836-BTA-51861 


21 


20990602 


ABHD2/PUN1 


8117/512222 


1 .39E-06 


0.056 
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Table 6. Cont. 





Trait 


Rank** 


SNP 


Chr. 


Position(bp) 


Nearest Gene/ 
Candidate Gene 


Distance(bp) 


Raw 

P_value' 


P_value 
Bonferroni 


UFA 


6 


Hapmap44836-BTA-51861 


21 


20990602 


ABHD2/PUN1 


8117/512222 


2.09E-06 


0.085 


UFA 


32 


Hapmap35708- 
SCAFFOLD316799_27843 


21 


25136856 


SH3GL3 


52052 


2.21 E-05 


0.896 


SFA/UFA 


39 


Hapmap33890- 
BES3_Contig41 8_1 1 54 


23 


402361 75 


ATXN1 


Within 


2.39E-05 


0.968 


SFA 


28 


ARS-BFGL-BAC-28144 


25 


2606575 


LOC788915 


21466 


1 .96E-05 


0.795 


SFA/UFA 


7 


ARS-BFGL-BAC-28144 


25 


2606575 


LOC788915 


21466 


2.89E-06 


0.117 


UFA 


5 


ARS-BFGL-BAC-28144 


25 


2606575 


LOC788915 


21466 


2.06E-06 


0.083 


SFA 


26 


BTA-61650-no-rs 


26 


41850719 


FGFR2 


Within 


1.91 E-05 


0.776 


SFA/UFA 


35 


BTA-61650-no-rs 


26 


41850719 


FGFR2 


Within 


1 .98E-05 


0.802 


SFA/UFA 


33 


BTB-01 926888 


27 


16398882 


TRIML2/ACSL7 


303782/2110549 


1 .96E-05 


0.797 


UFA 


21 


BTB-01 926888 


27 


16398882 


TRIML2/ACSL7 


303782/2110549 


1.58E-05 


0.641 


SFA 


18 


ARS-BFGL-NGS-1 06901 


29 


4437261 1 


SCYLl 


Within 


1 .22E-05 


0.495 


SFA/UFA 


30 


ARS-BFGL-NGS-1 06901 


29 


4437261 1 


SCYLl 


Within 


1.63 E-05 


0.663 



Note: see note to Table 2. 

doi:1 0.1 371/journal.pone.0096186.t006 



on BTA21 and Hapmap49848-BTA-106779 on BTA5, respec- 
tively. The SNP strongly associated with ClOiO (P= 8.54E-06), 
C12:0 (P=1.16E-07) and C14:0 (P= 6.01E-06), ARS-BFGL- 
NGS-39328, is 58,172 bp close to the fatty acid synthase [FASN] 
gene on BTA19, which is well-known to affect fat composition of 
dairy cattle and beef 

Long-chain saturated fatty acid traits (LCFA) 

A total of 126 significant SNPs for LCFA were detected mainly 
on BTAl, 2, 8, 10 and 17, including 105 for C18:0, 18 for C16:0, 
2 for C20:0 and one for C22:0. The top one significant SNP (ARS- 
BFGL-BAC-13788) was associated with C18:0 (P=9.17E-09) on 
BTA12. The strongest association of C22:0 (P=6.70E-07) was 
identified with the SNP (ARS-BFGL-NGS-1 09692) on BTAl. The 
SNP (ARS-BFGL-NGS-4939) associated with C16:0 (P=9.04E- 
06) on BTAl 4 is located within the diacylglycerol O-acyltransferase 
1[DGAT1) gene, the major gene with large effect on mUk fat in 
dairy cattle. 

Monounsaturated and polyunsaturated fatty acid traits 
(iVlUFA and PUFA) 

A total of 93 and 13 Significant SNPs for MUFA and PUFA 
were detected, respectively. Of them, 67, 6 and 20 SNPs were 
associated with CI 4:1, CI 6:1 and C18:ln9c, respectively. For 
C14:l, 29 out of 34 genome-wise significant SNPs on BTA26 were 
clustered within three regions: 6 fell in a 6.40 Mbp region 
(10.39-16.79 Mbp), 16 feU in a 1.75 Mbp region 
(20.36~22.11 Mbp) containing the stearqyl-CoA desaturase [SCDl) 
gene, and 7 SNPs fell in a 4.37 Mbp region (22.88-27. 25 Mbp). 
The SNP (BTB-00931481) on BTA26, 78,088 bp near to the 
SCDl gene, showed the strongest effect (P= 7.08E-13). Though no 
SNPs for CI 6:1 and C18:ln9c reached genome-wise level, the 
second suggestive significant SNP (ARS-BFGL-NGS-4939) for 
C18:ln9c (P=3.01E-06) is located within the DGATl gene on 
BTAl 4. 

As for PUFA, 6 and 7 significant SNPs were detected for 
C18:2n6c and CLA, respectively. The most significant SNP (ARS- 
BFGL-NGS-4939) associated with C18:2n6c (P=6.87E-08) is 
located within the DGATl gene on BTAl 4, while the most 



significant SNP (BTA-50525-no-rs) for CLA only reached 
suggestive level (P= 1.16E-05). 

Indices of fatty acid traits (CI 4 index, CI 6 index, CI 8 
index) 

For indices of C14, C16 and C18, totally 84, 10 and 93 
significant SNPs were detected, respectively. Forty-two SNPs 
associated with C14 index are located within a region of 
16.89 Mbp on BTA26, which included four small segments: 6 in 
a 2.40 Mbp segment (9.86- 12.26 Mbp), 7 in a 4.84 Mbp 
segment (14.15 — 18.99 Mbp), 18 in an 1.75 Mbp segment 
(20.36-22.11 Mbp) containing the SCDl gene, and 11 in a 
3.87 Mbp segment (22.88-26.75 Mbp). In addition, 56 common 
SNPs for C 14 index and CI 4:1, 4 common SNPs for C 18 index 
and C18:ln9c, 28 common SNPs for C 18 index and CI 8:0, and 4 
common SNPs for C16 index and C16:lwere identified. 

Sum of fatty acid traits (SFA, UFA, SFA/UFA) 

A total of 108 significant associations mainly on BTA5, 10 and 
20 with three sum of fatty acid traits were detected, which involved 
52 distinct SNPs. Of them, 22 SNPs were simultaneously 
associated with three traits and 12 were common for two traits. 
The 0.96 Mbp region (72.69-73.65 Mbp) on BTAIO was 
associated with the three traits, in which the SNP (ARS-BFGL- 
NGS-4783) showed the strongest association for SFA {P= 
6.07E-08), UFA (P=4.01E-07) and SFA/UFA (P= 2.64E-07), 
respectively. 

Discussions 

To our knowledge, this is one of the first GWA study for milk 
fatty acids with high density SNP Chip. In this study, we detected 
a total of 83 genome-wise and 314 suggestive significant SNPs for 
22 fatty acid traits. Among them, some SNPs are located within 
the QTL regions on BTA6, 14, 19 and 26 those have been 
reported by Stoop et at [13], Schenninkei al [14] and Morris et al 
[30] for bovine milk fat composition. Sixteen SNPs on BTAl 4, 5 
SNPs on BTAl 9, and 5 SNPs on BTA7 were consistent with the 
previous GWA study for fatty acid traits of dairy cattle [24]. 
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Table 7. Numbers of significant SNPs with genome-wise and 
suggestive significance for 18 milk fatty acid traits. 



Trait 


Genome-wise 
level 


Suggestive 
level 


Total 


CI 0:0 


1 


20 


21 


C12:0 


2 


20 


22 


CI 4:0 


1 


26 


27 


CI 4:1 


34 


33 


67 


CI 6:0 


0 


18 


18 


C16:l 


0 


6 


6 


CI 8:0 


13 


92 


105 


C18:ln9c 


0 


20 


20 


C18:2n6c 


3 


3 


6 


CLA 


0 


7 


7 


CI 4 index 


49 


35 


84 


C16 index 


0 


10 


1 0 


CI 8 index 


14 


79 


93 


SFA 


2 


28 


30 


UFA 


2 


37 


39 


SFA/UFA 


3 


36 


39 


C20:0 


0 


2 


2 


C22:0 


1 


0 


1 


Sum 


125 


472 


597 



Note: Associations with C8:0, C18:3n3, C18:3n6 and C20:5n3 only reached 
chromosome-wise significance or non-significant, so they were not listed. 
doi:l 0.1 371/journal.pone.00961 86.t007 



However, associations of BTA19 with CI 6:1 and CLA were not 
found in this study. This is probably due to a different dairy 
population was tested. Several SNPs were found to be located 
within and/or close to genes that are known to have functions 
related to the milk composition. In addition, 20 novel prospective 
candidate genes affecting milk fatty acid traits were identified. 

Chromosomes underlying novel promising candidate 
genes 

On BTAI, 23 SNPs associated with 9 fatty acids (CIO:0, C12:0, 
C14:0, CLA, C18:0, CIS index, C22:0, SFAand SFA/UFA) were 
detected. The SNP associated with SFA and SFA/UFA is 
38,733 bp away from the 3-hjdroxjacjl Coenzyme A dehydrogenase 
(EHHADH) gene. As a bi-functional enzyme, EHHADH is part of 
the classical peroxisomal fatty acid P-oxidation pathway, which is 
highly inducible via peroxisome proliferator-activated receptor a 
(PPARot) activation [31] and is essential for the production of 
medium-chain dicarboxyhc acids [32]. Four SNPs for CI8:0 and 
C18 index form an 0.40 Mbp region containing the 1-acylglycerol - 
3- phosphate O-acyltransferase 3 {AGPAT3) gene. AGPAT catalyzes 
the first step during de novo synthesis of triacylglycerol. AGPAT3 
is a member of the acyl transferase family [33] and plays a key role 
in de novo phospholipid biosynthetic due to its function of 
converting lysophosphatidic acid into phosphatidic acid [34]. 

On BTA2, 21 SNPs showed associations with 7 traits (CIO:0, 
C12:0, CI4:0, CI4:I, C18:0, CI8 index and UFA). Two SNPs 
associated with CI 8:0 and CIS index are 0.50 Mbp away from the 
signal transducer and activator of transcription 1 (STATl) gene, especially, 
one of them is the top 3 significant SNP for C 1 8 index. STATs are 
transcription factors known to importance to cytokine signaling. 



STAT 1 has a role in regulating the transcription of genes involved 
in milk protein synthesis and fat metabolism in Holstein [35] . The 
SNP associated with UFA is 0.14 Mbp and 19,295 bp away from 
the acyl-CoA synthetase long-chain family member 3 {ACSL3) gene and 
the monoacylglycerol O-acyltransferase 1 (MOGATl) gene, respectively. 
ACSL3 is an isozyme of the long chain fatty acids coenzyme A 
ligase family that convert free long chain fatty acids into fatty acyl- 
CoA esters and has a substrate preference for PUFA [36]. 
Depletion oiACSL3 by RNAi causes a significant reduction in fatty 
acids uptake, thereby plays a key role in lipid biosynthesis and fatty 
acids degradation [37]. MOGATl catalyzes the synthesis of 
diacylglycerols, the precursor of triacylglycerol and phospholipids 
[38], Two SNPs associated with CI8:0 and CI8 index are 
5.70 Mbp and 5.74 Mbp away from the fatty acid binding protein 3 
[FABP3) gene, respectively. FABP3 provides fatty acids for SCD, 
which is one of specific transporters for LCFA and one of the most 
abundant isoforms in bovine mammary tissue [39]. Eight 
contiguous SNPs associated with CI8:0 and CIS index are 
located within a chromosome region of 63.5S~98.16 Mbp that 
overlaps a reported QTL region (67.56~6S.25 Mbp) for CI4 
index, C 16 index, C18 index, SFA, MUFA, PUFA and SFA/UFA 
[40]. 

On BTA5, 17 SNPs showed association with 9 traits (C10:0, 
C14:0, CI6:0, C18:0, CI8:ln9c, CIS index, SFA, SFA/UFA and 
UFA). The top one significant SNP for CI4:0 is within the 
carboxypeptidase M [CPM] gene. Up-regulation of CPM in macro- 
phages (MAs) is associated with increased lipid uptake [41] and the 
highest expression of CPM was detected in human adipocyte cell 
[42]. The SNP associated with C18:In9c, SFA, UFA and SFA/ 
UFA is 0.77 Mbp away from the oxysterol binding protein-like 8 
[0SBPL8) gene which encodes a member of the oxysterol-binding 
protein (OSBP) family, a group of intracellular lipid receptors. 
0SBPL8 has the capacity to modulate lipid homeostasis and 
SREBP activity probably through an indirect mechanism [43] and 
is a negative regulator of sequestering of triglyceride [44] . The 
chromosome region of 8.43 Mbp (95.74~I04.17 Mbp) associated 
with C18:In9c, CIS index, UFA and SFA/UFA contains the 
oxidized low density lipoprotein (lectin-like) receptor 1 (OLRl) gene which 
can bind and degrade oxidized low-density lipoprotein [45,46] . 

On BTA9, 13 SNPs showed association with 9 traits (CIO:0, 
C12:0, CI4:0, CI6:I, CIS:0, CIS index, CIS:ln9c, CLA and 
UFA). The top one significant SNP for CIO:0 and C12:0 and the 
SNP for C18:0 are 0.40, O.Sl Mbp away from the 5 -hydroxy tryp- 
tamine (serotonin) receptor IB [HTRIB] gene, respectively. HTRIB is 
one of receptors for 5-hydroxytryptamine (serotonin). HTRIB 
gene knocked-out mice showed elevated aggression, higher food 
intake and impulsivity, indicating it possibly acts as a bridge 
between behavior and energy homeostasis [47]. Fatty acids, as 
energy signal, affect the activity of hypothalamic fat-sensitive 
neurons and impair nervous control of energy homeostasis [48] . 
HTRIB was also shown to affect milk production performance in 
Chinese Holstein [49]. 

On BTA20, 30 SNPs showed association with II traits (CI 2:0, 
C14:0, CI4:I, C18:0, CIS index, CIS:In9c, CIS:2n6c, CLA, 
SFA, SFA/UFA and UFA). The SNP associated with CI 8:0 is 
located within the growth hormone receptor (GHR) gene, the well- 
known major gene affecting milk fat trait [50]. Three SNPs 
associated with CIS:0 and CIS index are within an 1.46 Mbp 
region containing the 3-oxoacid CoA transferase 1 (OXCTl) gene 
which has a major function to utilize ketone bodies by mammary 
[39]. The SNP associated with SFA, UFA and SFA/UFA is 
located within the prolactin receptor (PRLR) gene which activates the 
STAT5A {Signal transducer and activator of transcription 5A) expression 
[51] and is associated with milk composition traits [52]. 
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On BTA21, 10 SNPs showed association with 8 traits (CI 0:0, 
C12:0, CI 6:0, C18 index, CLA, SFA, SFA/UFA and UFA). The 
0.15 Mbp region (9.37-9.52 Mbp) associated with C10:0 and 
C12:0 is 1.10 Mbp away from the insulin-like growth factor 1 receptor 
(IGFIR) gene. Furthermore, the top one significant SNP for CI 2:0 
is located within such region. IGFIR was found to affect milk 
composition traits [53]. Two SNPs associated with C10:0, SFA, 
SFA/UFA and UFA are close to the lipin 1 [PLINl) gene which 
plays a vital role in regulation on the expression of genes involved 
in milk fat synthesis [39] . 

On BTA26, 71 SNPs showed association with 8 traits (C10:0, 
C12:0, C14:l, C14 index, C16:0, C18:0, SFA and SFA/UFA). 
The nearest SNP is 32,576 bp close to the SCDl gene which 
encodes key enzyme responsible for the conversion of SFA to 
MUFA in mammalian adipocytes [54] and were shown to be 
associated with milk fetty acids [10,11,12,40,55]. The SNP 
associated with C10:0 and C14 index is located within the protein 
kinase, cGMP-dependent, type I [PRKGl) gene which is a key regulator 
of adipokine secretion and browning of white fat depots [56] and 
brown fat cell differentiation [57]. The SNP associated with C12:0 
is located within the multipk inositol-polyphosphate phosphatase 1 
{MLNPPl) gene. MlMPPl encodes multiple inositol polyphosphate 
phosphatase which converts 2, 3 bisphosphoglycerate (2,3-BPG) to 
2-phosphoglycerate [58]. As known, 2,3-BPG is a key substrate for 
the triglyceride (TG) synthesis. Two SNPs associated with CI 4:1 
and C14 index are 7,406 bp and 35,288 bp away from the lipase, 
family member J (LIPJ) gene and the lipase, family member K {LIPK) 
gene, respectively. The two genes belong to lipase family and take 
part in lipid catabolic process in human [59] and have an essential 
function in lipid metabolism of the most differentiated epidermal 
layers. The SNP for CI 4:1 and \he one for C 14:1 and C 14 index 
are 0.71 Mbp and 0.89 Mbp away from the enoyl coenzyme A 
hydratase, short chain, 1{ECHS1) gene, respectively. ECHSl takes 
part in fatty acid biosynthesis, elongation and metabolism and 
catalyzes the P-oxidation of fatty acid in human [60] . Two SNPs 
associated with C14:l and C14 index are located within the sorbin 
and SH3 domain containing 1 (SORBSl) gene. SORBS 1 is an 
important protein in the insulin-signaling pathway in the adipose 
depots of human [61] and has a positive regulation of lipid 
biosynthetic process [62]. The SNP associated with CI 4:1 and 
C14 index is located within the conserved helix-loop-helix ubiquitous 
kinase (CHUK) gene. CHUK takes part in mammary gland alveolus 
development, mammary gland epithelial cell proliferation [63] and 
Upogenesis through NF-kB activation pathway [64]. The SNP 
associated with CI 4:1 and C14 index is 1,586 bp away from the 
nuclear factor of kappa light polypeptide gene enhancer in B-cells 2 (NFKB2) 
gene which is essential for normal development of the mammary 
gland [65]. 

Chromosomes underlying known candidate genes 

Apart from the seven chromosomes as mentioned above, other 
chromosome regions harboring several significant SNPs within or 
near to known genes involved in fatty acids synthesis were 
identified. 

On BTA14, the SNP associated with 9 traits (C14:l, C16:0, 
CI8:ln9c, C18:2n6c, C14 index, C18 index, SFA, UFAand SFA/ 
UFA) is located within the DGATl gene, which has been 
confirmed to be the true QTL for milk fat composition in dairy 
catde [66]. On BTA19, three SNPs associated with C10:0, C12:0 
and C14:0 are 0.05, 1.71 and 1.42 Mbp away from FASN, ACACA 
(Acetyl-CoA carboxylase alpha) and SREBFl (Sterol regulatory element 
binding transcription factor 1), respectively. FASN is a multifunctional 
enzyme with a central role in the de novo lipogenesis in mammals 
[45,67]. ACACA catalyses biosynthesis of LCFA in mammalian 



cytosol [68]. SREBFl is a transcription factor that regulates the 
expression of the SCDl gene which is related to several genes of 
lipid metabolism [69]. On BTA6, the SNP associated with 

C 1 8index is located within the peroxisome prolferator-activated receptor 
gamma, coactivator 1 alpha (PPARGCIA) gene which is involved in the 
regulation of fatty acids transcription and mammary gland 
metabolism [70]. Six SNPs on BTA13 associated with C16:l, 
C16 index, C18:0 and C18 index are located 1.01-3.11 Mbp 
away from the acyl-CoA synthetase short-chain family member 2 (ACSS2) 
gene. ACSS2 provides activated acetate for de novo fatty acids 
synthesis [39]. Three SNPs within an 2.34 Mbp segment 
(54.27-56.61 Mbp) on BTA18 associated with C10:0, SFA, 
SFA/UFA and UFA harbors the Sphingosine kinase 2 {SPHK2) gene. 
As a lipid mediator with both intra- and extracellular functions, 
SPHK2 has diacylglycerol kinase activity and involves the 
sphingolipid synthesis [39] . Two SNPs on BTA23 associated with 
C16:0 is located 1.21 Mbp and 0.79 Mbp away from die 
butyrophilin, subfamily 1, member Al (BTNlAl) gene, which is essential 
for mUk hpid droplet formation [39], and the PRL gene, which 
impacts milk fat composition through STAT5A [45], respectively. 

The 0.03 Mbp region (16.39-16.42 Mbp) on BTA27 associat- 
ed with CI 2:0, SFA/UFA and UFA is 2.1 1 Mbp away from die 
acyl-CoA synthetase long-chain family member 1(ACSL1) gene. ACSLl 
has a vital role in fatty acids activation for milk TAG [39]. The 
SNP for C14 index is 1.12 Mbp away from the l-acylglycerol-3- 
phosphate 0-acyltransferase 6 (AGPAT6) gene, a novel lipid biosyn- 
thetic gene required for triacylglycerol production in mammary 
epithelium, i£ AGPAT6 was knocked out, lactating mice failed to 
synthesize milk fat [71]. The SNP on BTA29 associated with C18 
index is located 0.29 Mbp away from the fatty acid desaturase 
1{FADS1) gene which catalyzes the .synthesis of LCFA [72]. 

No significant SNPs were detected with C8:0, C18:3n3, 
C18:3n6 and C20:5n3, probably because these four traits have 
special population requirements. 

Conclusions 

The present genome-wide association study idimtificd 83 
genome-wide and 3 1 4 suggestive significant SNPs associated with 
18 milk fatty acid traits. Some of these SNPs were located within 
or near to previously reported genes and QTL regions, while some 
of the SNPs were novel. Consequentiy, 20 novel promising 
candidate genes were identified for C10:0, C12:0, C14:0, C14:l, 
C14 index, CI 8:0, C18:ln9c, C18 index, SFA, UFA and SFA/ 
UFA, such as HTRIB, CRM, PRKGl, MINPPl, LIPJ, LIPK, 
EHHADH, MOGATl, ECHSl, STATl, SORBSl, MFKB2, AGPAT3, 
CHUK, 0SBPL8, PRLR, IGFIR, ACSL3, GHR and OXCTl. Our 
findings are helpful for follow-up studies to fine-mapping to 
unravel causal mutations for milk fatty acid traits in dairy cattie. 
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Figure SI Manhattan plots for each studied milk fatty 
acids trait. BTAX is represented by BTA30), the first line 
represents genome-wise significant level (raw P<1.23E-06), and 
the second line represents suggestive significant level (raw P< 
2.46E-05). 
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